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FOREWORD 


The  research  was  performed  under  Program  Element  62601F,  Project  8809, 

Task  04. 

Inclusive  dates  of  research  were  1  September  1971  through  1  October  1973. 
The  report  was  submitted  25  October  1973  by  the  Air  Force  Weapons  Laboratory 
Project  Officer,  Major  Arthur  R.  Geldbach  (SA3). 

The  advanced  computer  model  QUANTO  has  been  developed  within  the  Air  Force 
Weapons  Laboratory  to  study  various  scenarios  involving  sea-launched  ballistic 
missile  attacks  cn  bomber  air  bases.  The  QUANTO  model  has  been  reviewed  by 
interested  Air  Staff  agencies,  the  Air  Force  Systems  Command,  and  the  Strategic 
Air  Command,  and  is  considered  appropriate  for  use  in  activities  relating  to 
bomber  force  prelaunch  survival.  However,  prudence  should  be  exercised  in  its 
use,  because  of  its  sensitivity  and  the  dynamic  nature  of  the  problem. 

The  basic  model  was  devel :oed  by  Major  Richard  Conway.  A  large  portion  of 
the  debugging  and  exercising  o:  the  code  was  done  by  Mr.  Eugene  Omoda  and  Mr. 
William  Peay.  The  assistance  of  Mr.  Harry  Murphy  in  uti’izing  the  operating 
system  and  remote  terminal,  and  that  of  Mr.  A1  Sharp  in  incorporating  the 
thermal  and  overpressure  routines  into  QUANTO  were  also  invaluable  in  the 
development  process. 

This  technical  report  has  been  reviewed  and  is  approved. 


ARTHUR  R.  GELDBACH 
Major,  USAF 
Project  Officer 


.  JttRGf  H.  DIMON^R^*^" 
p  Colonel ,  USAF 

Chief,  Battle  Environments  Branch 


CLLc 


CHARLES  C.  HYRE,  JR 
Colonel,  USAF 
Chief,  Analysis  DVFision 
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SECTION  I 
INTRODUCTION 

The  theory  of  the  allocation  of  the  sea-launched  ballistic  missiles  (SLBM) 
against  a  force  of  aircraft  flushing  from  their  respective  airbases  and  the 
defensive  reactions  to  given  threat  levels  is  discussed.  The  analysis,  which 
led  to  the  models  discussed  later,  has  culminated  in  a  computer  program  called 
QUANTO.  The  model  used  in  QUANTO  has  as  its  inputs  latitude  and  longitude 
coordinates  of  target  and  submarine  locations,  aircraft  beddowns,  aircraft  and 
missile  flight  parameters,  and  aircraft  vulnerability  levels.  Consequently, 
the  code  is  useful  for  studying  the  effects  of  variations  in  a  number  of 
parameters. 

QUANTO  analyzes  three  types  of  problems  important  to  strategic  planners: 

Case  I:  Given  specific  locations  (Aa)  for  a  fixed  number  of  attacking 
submarines  and  t  specific  beddown  for  aircraft  at  locations  T^ ,  QUANTO  can 
compute  where  the  assigned  missiles  from  A^  should  go. 

a 

Case  II:  Given  specific  beddown  for  aircraft  at  locations  T.j,  QUANTO 
can  optimize  the  locations  for  the  submarines  among  a  set  of  candidate  locations 


I  • 


Case  III:  Given  specific  subma.  ne  locations  A  .  QUANTO  can  optimize 

a 

the  beddown  of  aircraft  at  T^. 

Lagrange  multipliers  are  used  in  the  optimization  procedures  of  QUANTO.  A 
br  ef  review  of  this  technique  is  presented  in  appendix  I  and  is  intended  to 
acquaint  the  reader  with  the  basic  mathematics  involved. 

The  QUANTO  code  has  been  developed  within  the  Air  Force  Weapons  Laboratory 
(AFKL).  It  was  intended  originally  as  a  vehicle  for  increasing  the  understanding 
of  the  operation  of  a  computer  program  called  COG,  which  dealt  only  with  Case  I 
(as  uf  May  1971),  that  was  written  by  the  Lambda  Corporation  (ref.  1).  Compared 
to  other  codes,  QUANTO  permits  a  more  detailed  and  accurate  analysis,  because 
weapons  and  their  detonations  arc  handled  individually,  rather  than  as  members  of 
fixed  weapon  patterns.  Studies  show  that  QUANTO  produces  a  considerably  better 
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allocation  than  does  COG.  Further  investigation,  substantiated  by  simulation 
of  the  attack  through  the  use  of  another  AFWL  code,  supports  the  assumptions  and 
models  used  in  QUANTO.  Hence,  QUANTO  provides  a  means  for  comparing  and  evalu¬ 
ating  the  effectiveness  of  other  weapon  allocation  codes.  More  importantly, 
QUANTO  provides  a  framework  for  modification  and  extension  in  further  studies 
of  total  bomber/tanker  force  survivability. 
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SECTION  II 

BASIC  WEAPON  ALLOCATION  PROBLEM 

In  this  study,  the  attacking  force  of  submarines  (figure  1)  is  distributed 
among  points  A  ,  a  =  ],  2,  p,  where  submarines  at  A3  each  carry  s,  SL8M 

«  a  a 

weapons.  At  the  same  time,  suppose  that  the  targets,  T. ,  i  =  1 ,  2,  . . . ,  M, 
have  values  V..  If  one  were  to  visualize  this  engagement  as  in  figure  1,  it 
becomes  apparent  that  many  strategies  are  open  to  the  a^tacher  and  defender. 

For  example,  the  attacking  force  could  put  all  missiles  on  target  Tj.  On  the 
other  hand,  the  missiles  could  be  distributed  among  all  targets.  As  for  the 
defender,  he  could  place  his  bombers  and  tankers  throughout  Lha  target  areas 
evenly  or  perhaps  all  on  the  same  base.  The  multiplicity  of  possibilities 
increases  with  each  new  missile  or  aircraft,  making  hand  calculation:-)  impracti¬ 
cal.  The  approach  taken  to  solve  this  problem  is  to  use  the  method  of  Lagrange- 
Multipliers  to  produce  a  near-optimal  allocation  of  SLBMs  to  targets  consisting 
of  escaping  aircraft.  To  construct  the  objective  function  which  describes  the 
expected  value  killed,  one  must  first  develop  the  survival  probability  S^.  for 
each  weapon.  This  figure  is  the  probability  that  a  single  aircraft,  given  an 
escaping  time-dependent  pattern  of  aircraft  taking  off  from  the  airfield, 
survives  one  incoming  SLBM.  The  probability  of  kill  is  then 

■V-i-s,  <’> 

Suppose  now  that  n  weapons  are  delivered  to  a  target,  and  the  survival  probabil¬ 
ities  $.  of  the  target  from  each  weapon,  j  =  1,  2,  ...,  n,  are  independent;  then 

J 

the  probability  of  destroying  the  target  is 

Pk  =  1  -  S  S.  (2) 

k  j=l  0 

For  a  system  of  M  targets,  each  having  a  value  V..,  the  expected  return  from 
delivery  of  all  weapons  is 
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where  S^  is  the  probability  of  survival  of  target  i  from  weapon  j.  ihe  product 
in  expression  (3)  for  target  i  includes  only  those  survivabilities  corresponding  to 
the  weapons  which  are  aimed  at  target  i.  In  practice,  there  are  weapon  "groups," 
where  the  weapons  in  each  group  are  so  nearly  identical  in  characteristics  and 
location  that  no  distinction  between  them  is  required  for  purposes  of  allocation. 
Hence,  in  practice,  S. .  is  raised  to  the  n--  power,  where  n- .  is  the  number  of 

•  J  •  V*  *  J 

weapors  from  weapon  group  j  that  are  targeted  against  target  i.  It  is  easy  to 
see  that  the  product 

n. . 

n  $..  13 

V  10 

J 

is  not  changed  if  the  weapon  groups  j  are  included  for  which  n,-  •  =  0.  Conse- 

•  J 

quently,  if  L  is  the  total  number  of  weapon  groups,  the  expected  aircraft  kill 
may  be  written 


£  v,  i  -  n  s,, 
i=l  1  L  j=l  10 


Table  I  clarifies  the  submarine  input  parameters  used  in  QUANTO.  The  table 
has  six  columns,  the  first  of  which  is  submarine  locations,  given  to  QUANTO  in 
terms  of  latitude  and  longitude  coordinates,  surrounding  a  given  target  country. 
The  second  and  fifth  columns  contain  the  same  information  and  are  presented 
separately  to  emphasize  the  fact  that  the  number  of  submarines  and  number  of 
missiles  per  group  are  the  same  since  all  the  submarines  at  a  given  location  are 
assumed  to  fire  a  missile  apiece  at  the  same  time.  Note  here  that  zero  sub¬ 
marines  are  allowed  at  a  given  location.  In  column  four  the  numbering  1  to  4  is 
applied  to  two  types  of  weapons,  each  of  which  is  restricted  to  either  the 
Atlantic  or  Pacific  Ocean.  Numbers  1  and  2  may  identify  weapon  types  1  and  2  in 
the  Pacific,  whereas  3  and  4  may  represent  weapon  types  1  and  2  in  the  Atlantic. 
Since  submarines  may  be  shifted  only  among  locations  which  have  like  missile 
type  identifiers  in  QUANTO,  such  a  numbering  system  prevents  submarines  from 
relocating  to  a  different  ocean.  The  last  column  is  the  numbering  given  to  the 
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Table  I 

SUBMARINE  INPUT  DATA  TO  QUANTO  (EXAMPLE) 


Location 

Number 

Submarines 

Number 

Missiles/ 

Sub 

(L  =  42) 

Type 

Missile 

Number 

Missiles/ 

Group 

Cj> 

Weapon 

Groups 

(j) 

1 

1 

8 

1 

1 

1  -  3 

2 

n 

C. 

6 

2 

2 

9  -  14 

3 

0 

8 

1 

0 

15  -  22 

4 

4 

6 

2 

4 

23  -  28 

5 

1 

8 

3 

1 

29  -  36 

6 

3 

6 

4 

3 

37  -  42 

weapon  groups.  Note  here  that  the  numbering  in  row  2  goes  from  9  through  14. 
Each  weapon  group  here  consists  of  two  missiles  in  the  same  salvo  since  there 
are  two  submarines  at  this  location.  Missiles  may  be  placed  in  the  same  group 
if  they  have  identical  trajectories  and  are  launched  at  the  same  time  from  the 
same  point.  Also,  row  3  allows  for  a  set  of  weapon  groups  even  though  no  sub¬ 
marines  are  initially  placed  at  submarine  location  3  (although  there  may  be 
subsequently,  if  the  submarine-placement  optimizer  of  Case  II  is  exercised). 

The  basic  allocation  problem  is  to  maximize  the  expected  kill  given  by 
expression  (4.  by  sending  the  missiles  to  the  proper  targets.  Since  the  alloca¬ 
tion  of  missiles  to  targets  is  expressed  by  the  values  n. .,  the  problem  is  to 

*  vl 

find  the  integer  values  n*.  which  result  in  the  greatest  kill  while  satisfying 

*  J 

constraints  on  the  number  of  weapons  available  in  each  group. 
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SECTION  III 

USE  OF  LAGRANGE  MULTIPLIERS  IN  THE  QUANTO  COMPUTER  CODE 


1.  FORMULATION  FOR  SOLUTION 

The  weapon  allocation  problem  is  one  of  determining  the  optimal  allocations 
n-  of  weapons  to  targets  to  maximize  the  expected  kill  value 

*  V 


v  M  r  L  n  ■ 

K)  *  £  i1  -  °  ,Jj 


subject  to  the  stockpile  constraints 


£  "u =  V 


=  12  L 

«-  }  t  »  «  )  L 


where  N^  is  the  number  of  weapons  in  weapon  group  j.  Fractional  allocations  in 
n^  are  permitted  in  the  solution  of  this  ^oblem,  but  each  n^  must  satisfy 

(7> 

As  described  in  appendix  I,  this  problem  gives  rise  to  the  Lagrangian  function 

h(n'j’Xj)  =  f(nij)  +  £  [j?,  (nij)  ■  No]  (8) 

In  seeking  the  extremum  of  the  Lagrangian  function  h  n..,A.  ,  the  values  of 

(\  J  J 

n. are  sought  which  satisfy  the  following  necessary  conditions  for  a 
solution  using  this  Lagrangian  Multiplier  technique 


ah  _ 


=  -vk  (*»  skt)  n  +  X,  •  0 

J  ® 


k  =  1,  2,  M;  z  =  1,  2,  L 
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If  variable  x.  (dependent  on  n,,.-,  j  *  1,  2 .  L)  are  defined  as 


xk£  H  "Vk 


(£n  5u)  £ 


,  L  n. 


the  system  (equation  (9))  of  (M  X  L)  equations  becomes 

Xk£,  =  "V  k  =  U  2,  ....  M;  £  =  1,  2,  ....  L 
Now  fix  £  and  consider  the  subsystem  of  M  equations 


xu  =  “X£ 


X2i  ~xz 


XM£  ~  "X£ 

A  word  is  in  order  concerning  notation.  In  equations  (11),  x£  is  one  of  the 
unknown  Lagrange  multipliers.  The  variables  xk£  (k  =  1 ,  2,  . . . ,  M)  are  computable 
if  one  has  the  values  of  nkj  (j  =  1,  2,  ...,  L).  The  technique  used  for  finding 
the  values  of  x£  and  nkj  (k  =  1,  2,  ...»  M;  j  =  1,  2,  ...»  L)  which  satisfy  the 
system  (equations  (11))  of  M  equations  takes  advantage  of  the  fact  that  all  the 
xk£  should  equal  the  same  quantity,  namely  -x£.  The  method  chooses  values  of 
nkj  iteratively,  subject  to  the  constraints,  so  that  the  values  of  xk£  (k  =  1, 

2,  ...,  M)  approach  a  single  value,  namely  -x£. 

2.  ITERATIVE  PROCEDURE 

An  initial  allocation  of  weapons  to  targets  n. .  is  input  to  QUANTO,  and  the 
variables  A^ .  are  computed.  Suppose  for  a  given  weapon  group  £  that  xk£  <  xm£ 
and  nk£  >.  0.0001.  Then  by  moving  an  appropriate  number  An  of  weapons  in  group  £ 
from  target  k  to  target  m,  xk£  and  xm£  may  be  made  more  nearly  equal .  Note 
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that  nki  must  be  initially  positive  or  there  would  be  no  weapons  to  shift.  In 
fact,  were  it  not  for  the  restriction  that  n.  may  not  be  reduced  to  a  negative 

1\A* 

amount  (i.e.,  An  <.  nk£),  xki?  and  xm£  could  be  made  equal  in  all  cases.  The 
value  of  An  which  would  make  the  new  values  of  xk£  and  xm£,  say  xk£  and  xm£, 
equal  is  the  value  of  An  which  satisfies 

S 

S. 

m  i) 

This  may  be  written  as 


~  Vm  (£n  Smt) 


.  nkl 

<-  nk2 

-  nk£“An 

1 

C 

*kl 

Sk£ 

... 

n  i 
ml 

n  « 

r 

nmC+An 
o  m£ 

1 

c 

’ml 

bm2 

m£ 

•••  ^mL 

‘kL 


mL  _ 


=  X. 


m£ 


X  =  X  S  -An  =  X  S  +An  -  X 
k£  ks.  k£  m£  m£  m£ 


Therefore, 


£n 


An  = 


_k£_ 

\n£ 


fn  (S 


k£ 


v) 


02) 


Since  An  is  not  permitted  to  be  so  large  that  (nk£  -  An)  becomes  negative,  the 
actual  number  of  weapons  shifted  is 


An  =  mm 


*liwv) 

I""’  *  (Sk*  W 


(13) 


This  shift  of  weapons  gives  rise  to  a  new  j..  and  new  X--.  Repeated  shifts 

(\  *  J  *  J 

Ak£*  Am£/*  £or  eac^  weaP°n  9r0UP  l>  t0  &e  eRual 
(for  those  targets  k  and  m  for  which  weapons  from  group  £  end  up  being 
allocated). 

Although  the  restriction  An  <n^  makes  it  impossible  to  force  the  equality 
of  every  pair  ^xk£,  xm£j,  the  preceding  choice  of  An  does  result  in  the  greatest 
increase  in  the  objective  function  which  can  result  from  such  a  shift  of  weapons 
in  group  £  from  target  k  to  target  m.  To  see  this  consider  the  function 
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r(in)  =  Vk  [.-Sk)"kl  •  Sk/“  ...  ...  SkL"kL] 

*  vm  ['V"”  ■  VV  -  sm/mt+4"  -  smLV]  <14> 

which  represents  the  kill  contribution  to  the  objective  function  f(n-jj) 
(equation  (5))  from  targets  k  and  m  after  An  weapons  are  moved  from  target  k  to 
target  m.  The  best  choice  of  An  is  where  r(An)  achieves  its  maximum  within  the 
interval  0  <.  An  <.  n^.  The  unrestricted  maximum  of  r(An)  occurs  where 

d  r(An)  _  n 


i  .e. , 


r'(“)  ■  V  -k. J"  ■  0 


An*  = 


1(xki,/xme) 

-n  (Sk£Sm£) 


If  this  value  is  greater  than  the  constrained  maximum  of  r(An)  occurs  at 
An  =  .  This  follows  from 

r' "» =  v  -  v  > 0 


r'  •  (*n>  =  xm  (en  Sk,)  sj 


(in  smt) 


s.  ~&n  <  o 

k£  “ 


for  all  An  in  the  range  0  £  An  <.  nk£  (since  0  <  S.^  <  1,  £n  <  0  and  >.  0 
for  all  i,  j).  Thus,  r(An)  appears  as  in  figure  2  or  figure  3.  Note  that  the 
curvature  is  always  downward  and  that  the  maximum  occurs  at  the  point  An*.  If 
the  situation  of  figure  3  occurs,  it  is  impossible  to  choose  An  -  An*  to  force 
xk£  and  xm£  to  be  equal.  Consequently,  equations  (11)  will  not  be  satisfied. 
However,  the  optimal  value  of  f(n-jj),  where  the  n..  are  constrained  by  equation 


I 


j» 

> 

4 
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Figure  2,  Constrained  Maximum 
4n*  *■  nk2 


Figure  3.  Constrained  Maximum 
4n*  »  nki 


(6)  and  equation  (7),  is  still  found  due  to  the  preceding  comments  concerning 
r(an).  Thus,  the  X- - ‘s  merely  serve  as  a  means  of  adjusting  the  n-.'s  to 

*  J  *  J 

approach  optimality.  The  optimal  n  - .  are,  of  course,  nonintegral  and,  therefore, 

*  J 

not  physically  possible.  Consequently,  the  optimal  nonintegral  allocation  is 
integer! zed  to  give  an  integral  allocation  which  satisfies  the  constraints. 

This  integerizacion  is  performed  for  each  weapon  group  j  by  rounding  those  n-.'s 

'  J 

with  the  largest  fractional  parts  up  and  rounding  the  remaining  n-^-'s  down.  Of 
course,  the  constraints 


1 

s 


(in  which  the  N-'s  are  integers)  are  satisfied  by  the  real  n.  . 's  before  integeri- 

J  '  g 

zation  and  must  be  satisfied  by  the  integer  n-.'s  also.  Therefore,  the  number  of 

*  J 

n-.'s  rounded  up  is  determined  so  that,  for  each  weapon  group  j,  the  sum  of  those 

*  J 

n-.'s  rounded  uo  and  those  rounded  down  equals  N-.  In  practice,  the  expected 

*  v  J 

kill  resulting  from  this  integeri zed  allocation  is  not  significantly  different 
from  the  expected  kill  computed  from  the  nonintegral  allocation,  since  the  dif¬ 
ference  in  kills  is  usually  only  a  fraction  of  an  aircraft.  Integeri zation  of 
the  optimal  nonintegral  allocation  need  not  produce  the  optimal  integral  alloca¬ 
tion,  but  it  does  produce  at  least  a  near-optimal  integral  allocation,  with  the 
difference  in  kills  being  the  upper  bound  of  how  far  from  optimal  the  kill  of 
the  integeri zed  allocation  could  be. 

It  has  been  indicated  hew  n^  and  nm£  may  be  adjusted  to  increase  the 
expected  kill  value  when  X}.£  <  x^  and  nk£  z.  0.0001  for  some  weapon  group  fc. 


* 
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In  practice,  j  tolerance  level,  e,  is  set  in  QUANTO,  so  that  convergence  is 
said  to  occur  when  xm£  -  e  <.  xk£  <.  xm£  for  all  n^  >_  0.0001  for  all  values  of  t, 
where  Xm£  •=  max  jA^},  (i  =  1,  2,  ....  M). 

Specifically,  the  weapon  group  t,  upon  which  each  allocation  adjustment  is 
based,  is  selected  in  a  cyclical  manner.  The  first  allocation  adjustment  is 
made  within  weapon  group  one  (t  =  1)  if  x^  <  xm£  -  e  where 


xk£  =  min  |xi£  such  that  ni£  >_  0.000lj,  i  =  1,  2,  ...»  M 


06) 


and 


xm„ 

mu 


=  max 
i 


1  =  1,  2, 


•  •  •  9 


M 


(17) 


If  this  situation  does  not  exist  for  s,  =  1,  successive  weapon  groups  are 
inspected  in  sequential  order  until  one  is  found  in  which  the  highest  X^£  exceeds 
the  lowest  xi£  with  a  corresponding  positive  allocation  (n.j£  >.0.0001)  by  more 
than  the  tolerance  e.  Successive  allocation  adjustments  are  accomplished  in  a 
repeating  cycle  through  the  values  of  i  {i.e.,  1,  2,  ...»  L,  1,  2,  ...,  L,  1,  2, 
...).  Convergence  occurs  when  all  weapon  groups  are  inspected  without  finding 
one  which  initiates  an  allocation  adjustment.  In  practice,  the  A- -  matrix  if 

*  J 

first  converged  to  a  tolerance  of  e  =  0.1,  then  e  =  0.01,  then  e  =  0.001  and 
so  forth,  with  r.he  final  tolerance  under  the  control  of  the  user.  This  process 
results  in  a  faster  overall  convergence  to  the  final  tolerance  level.  An 
additional  cutoff  of  the  convergence  occurs  if  a  given  number  (sp'.cified  by  the 
user,  say  100)  of  allocation  adjustments  are  performed  without  increasing  the 
kill  value  above  some  user-selected  amount  (say  0.01).  The  user  may  also 
simply  specify  a  maximum  number  of  allocation  adjustments  to  be  made. 


The  iterative  procedure  is  illustrated  in  the  flow  chart  of  figure  4,  with 
several  additional  details  appearing  in  the  figure.  When  xk£  <  xm£  -  e  is  found 
for  some  weapon  group  t.  An  must  be  computed.  When  only  one  type  of  aircraft  is 
considered  in  the  model.  An  is  computed  according  to  equation  (12).  However, 
when  more  than  one  type  of  aircraft  is  considered.  An  must  be  computed  by  a 
Newton  iterative  procedure,  described  in  the  mixed  force  allocation  problem 
discussed  in  section  VI. 
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The  amount  of  output  may  be  controlled  to  some  extent  by  a  control  variable 
set  by  the  user.  Output  on  each  iteration  may  be  suppressed  and  only  a  limited 
output  obtained  after  a  "subgroup"  of  iterations.  If  the  expected  kill  has  not 
increased  significantly  for  the  subgroup  of  iterations,  the  procedure  is  termin¬ 
ated.  Two  other  conditions  may  terminate  the  procedure,  as  shown  in  figure  4 
and  previously  described. 

3.  OPTIMIZATION  OF  SUBMARINE  LOCATIONS  IN  QUANTO 

In  the  weapon  allocation  problem,  the  Lagrange  multiplier  A.  represents  the 

J 

shadow  value  associated  with  weapon  group  j.  In  notation 


Hr  *  -»j  0  (18> 

3 

where  h(n-jj,  Aj)  is  the  Lagrangian  function.  Therefore,  increasing  N.  has  the 
instantaneous  effect  of  permitting  an  increase  in  h  (and,  therefore,  f)  at  the 
rate  of  -A.  units  of  f  per  unit  of  N-.  Thus,  one  can  get  some  feel  for  the 

J  J 

value  of  an  additional  weapon  in  group  j  ob  observing  the  magnitude  of  A.. 

J 

A  heuristic  rule  has  been  used  in  QUANTO  to  relocate  submarines  among  the 
input  candidate  submarine  locations  so  as  to  improve  the  expected  kill  value. 
The  value  of  a  submarine  at  a  given  location  bears  some  relation  to  the  magni¬ 
tudes  of  Akj  for  those  weapon  groups  j  corresponding  to  salvos  from  a  submarine 
at  that  location.  A  submarine  location  input  to  QUANTO  is  characterized  not 
only  by  its  geographical  coordinates  but  also  by  the  type  of  submarine  (and 
its  number  of  salvos,  s)  which  can  be  located  there. 

For  each  submarine,  the  quantity 

is  calculated  where  the  sum  is  over  the  s  values  of  j  corresponding  to  that 
iubmarine's  salvos,  and  each 

xkj  =  ml*n  {x1j  suc^  ^at  nij  —  0*0001  j,  i  =  1,  2,  ...,  M 

for  each  weapon  group  j  at  that  submarine's  current  location.  Similarly,  for 
every  other  location  at  which  that  submarine  can  operate,  the  quantity 
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is  calculated,  where 


Amj  max  i  -  1,  2,  . • » ,  H 

♦  I 

with  a  view  to  possibly  moving  the  submarine  to  a  location  where  it  can  be 
expected  to  kill  more  value.  When  there  are  several  types  of  submarines,  the 
type  of  submarine  moved  is  the  one  having  the  largest  average  difference 


V  xmj  -  £  xkjj>  xkj  such  that  nkj  ^°-0001 

Within  this  submarine  type,  the  submarine  relocated  is  the  one  corresponding 
to  the  lowest  quantity 


and  it  is  placed  in  the  location  having  the  highest  value  of 


£  x  . 
j  mj 


Relocation  of  a  submarine  in  QUANTO  is  accomplished  by  moving  the  s  missiles 
on  that  submarine  to  another  location.  Consequently,  n. .  is  increased  by  one 

*  vl 

missile  in  the  s  weapon  groups  j  corresponding  to  the  submarine  location  to  which 
that  submarine  is  moved  and  those  additional  SLBMs  are  assigned  to  targets  i 
having  the  largest  value  X- •  (for  each  j).  Similarly,  for  the  s  weapon  groups  j 

•  J 

from  which  a  missile  is  removed,  n..  is  reduced  for  the  targets  i  corresponding 

*  0 

to  the  lowest  x  —  's  until  a  total  of  one  missile  is  removed  from  each  weapon 
group.  In  this  way,  a  rational  guess  is  made  at  where  the  missiles  from  the 
relocated  submarine  should  go  in  order  to  obtain  an  initial  allocation  prior  to 
re-entering  the  laydown  optimization  procedure. 

This  submarine  relocation  process  in  no  way  guarantees  an  increase  in  value 
killed.  This  is  because  the  relocation  is  accomplished  by  moving  an  integral 
number  of  missiles,  not  a  An  computed  to  maximize  kill.  Also,  s  missiles  (not 
just  one)  are  moved  before  the  x . •  are  recomputed  from  the  new  n. .'s. 

I J  *  %) 
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3F*  |* 

I  I: 


Although  the  heuristic  rule  does  not  always  increase  the  kill,  experience 
with  the  procedure  reveals  that  the  kill  usually  increases  with  every  more  until 
a  decrease  occurs,  after  which  the  kill  varies  with  additional  moves  without 
significant  gains  or  losses.  Consequently,  submarine  moving  is  terminated  in 
QUANTO  after  the  first  move  which  results  in  a  decreased  kill. 

4.  BEDDOWN  OPTIMIZATION 

A  heuristic  routine  for  shifting  aircraft  from  base  to  base  has  been  supplied 
in  the  QUANTO  weapon  allocation  code  with  the  intent  of  determining  better  air¬ 
craft  beddowns  for  a  given  positioning  of  submarines.  The  procedure  shifts  air¬ 
craft  from  the  base  having  aircraft  value  greater  than  0.0001  with  the  lowest 
survivability  product 

L  nii 

n  S.. 
j=l  1J 

to  the  base  with  the  highest  survivability  product,  provided  the  losing  base 
starts  with  an  aircraft  value  greater  than  0.0001.  Otherwise,  the  survivability 
products  are  inspected  in  ascending  order  until  a  corresponding  aircraft  value 
greater  than  0.0001  is  found,  and  the  corresponding  base  is  selected  as  the 
losing  base.  The  amount  of  value  shifted  is  the  nonintegral  product 


aV  =  V, 


L  n  . 
n  S  .  mj 
j-1 


n  Ski 
j=l  kJ 


where  bases  m  and  k  are  those  having  the  highest  and  lowest  survivability  prod¬ 
ucts,  respectively,  where  Vj,  >  0.0001.  If  <  0.1,  then  all  the  value  is 
moved  from  base  k  to  base  m  regardless  of  the  aV  computed. 

A  with  the  heuristic  submarine  relocation  routine,  each  shift  of  beddown 
value  In  accordance  with  the  above  formula  does  not  guarantee  a  decrease  in  the 
overall  expected  kill  value,  although  the  general  trend  is  toward  a  lowering  of 
the  kill.  Occasionally,  an  overall  kill  increase  may  occur  as  the  result  of 
individual  shifts  considerably  before  the  process  has  exhausted  the  gains  to  be 
made  in  aircraft  surviving.  The  shifting  of  value  terminates  if  the  value  aV 
to  be  moved  is  less  than  0.05  (specified  by  a  program  statement),  at  which  point 
the  survivability  products  have  essentially  converged  and  the  beddown  is  not 
changing  significantly. 
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Shifts  of  aircraft  do  not  cause  a  recomputation  of  the  s-rvivabilities  S.. 
in  QUANTO,  since  only  rarely  does  che  computed  survivability  depend  upon  the 
number  of  aircraft  present  at  the  target.  The  methods  of  computing  the  surviva¬ 
bilities  S..  are  described  in  the  next  section  of  this  report. 

*  J 

If  both  beddown  optimization  and  optimization  of  submarine  locations  are 
requested  by  the  user  of  QUANTO,  the  beddown  optimization  is  performed  last.  Of 
course,  if  the  user  wants  the  submarines  to  have  the  last  move,  he  may  request 
beddown  optimization  only,  and  in  a  subsequent  run,  input  the  optimal  beddown 
and  request  submarine  optimization. 

The  beddown  optimization  procedure  shifts  nonintegrai  numbers  of  aircraft, 
and  thus  results  in  a  beddown  which  has  fractional  numbers  of  aircraft  at  the 
various  bases.  After  the  termination  of  aircraft  moves,  the  beddown  is  inte- 
gerized  along  with  the  missile  laydown  and  the  results  are  output.  Integeriza- 
tion  of  the  beddown  has  the  effect  of  increasing  the  kill  by  a  negligible  amount 
over  the  expectea  kill  computed  on  the  basis  of  nonintegral  beddown. 

The  following  discussion,  in  the  form  of  a  critiqued  proof,  is  presented  as 
a  partial  justification  for  the  heuristic  rule  for  improving  the  beddown. 

Suppose  (Vj,  V2,  . ..,  Vk,  ...,  V  ,  ...,  VM}  represents  the  values  of  the  aircraft 
bedded  down  on  the  M  bases  for  which  the  optimal  missile  laydown  is  [n-j  j]  and 


l 

n  s 

j=i 


kj 


L 

<  n 
j=i 


(20) 
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Next  suppose  the  beddown  is  changed  by  subtracting  some  small  value  e  >  0,  from 
Vk  and  adding  e  to  V’m.  Thus,  (V1#  V2,  ...,  Vk  -  e,  ...,  Vm  f  s,  VM) 

represents  the  new  beddown,  and  the  new  optimal  laydown  [ft^j  could  be  determined. 

The  new  value  surviving  in  the  new  beddown  is  then  '  j 

M  l  n,,  L  n  .  L  n.. 

New  Surviving  Value  =  ^  V .  n  S..  +  e  n  S  .  m3  -  n  S,  .  J  (21) 

1=1  1  j=l  ,J  "W  j=i  kj 

Now  if  e  is  sufficiently  small,  it  is  reasonable  to  exiect  that  [n^]  is  close 
to  [n-j-j],  so  that  j 


M  L  ft.-,-  M  L  n 

£  v.  n  s. .  J  >  £  V,  n  s..  13  =  Old  Surviving  Value  (23) 
1=1  1  j=1  1J  1=1  1  j=l  13 


since  [nij]  is  the  optimal  laydown  for  the  old  beddown  (Vj)  and  therefore 
minimizes  the  survivors.  Consequently,  a  shift  of  value  {sufficiently  small) 
from  base  k  to  base  m,  when 


L 

n 

j=i 


n  . 
c  rnj 

mj 


L 


n 

=1 


result.'  in  a  reduced  expected  kill.  Note  that  this  proof  does  not  indicate  the 
best  amount  of  value  e  to  shift,  but  merely  that  value  should  be  shifted  to 
bases  having  high  survivability  products  from  those  with  low  products. 

5.  LETHAL  AREA  DETERMINATION 

The  determination  of  lethal  area  (i.e.,  the  region  within  which  aircraft 
are  destroyed)  resulting  from  a  nuclear  weapon  detonation  is  an  integral  part  of 
tha  QUANTO  code.  The  lethal  areas  are  required  in  the  computations  of  the 
survivabilities.,  S..,  of  aircraft  flushing  each  target  area.  This  subsection 

•  J 

discusses  the  assumptions,  assertions,  models,  and  methods  used  in  the  lethal 
area  determination. 

The  nuclear  environment  created  by  the  detonation  of  nuclear  weapons  is 
discussed  in  AFSCM  500-1.  This  section  is  concerned  with  the  fireball  effect 
(thermal)  and  the  blast  (overpressure)  effect  which  are  considered  to  be  the 
only  two  structural  kill  mechanisms  which  can  destroy  an  aircraft  for  this  model. 
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A  superheated  region,  the  fireball,  cools  while  expanding  and  radiates  thermal 
energy  (heat).  At  the  expanding  edge  of  the  fireball,  tremendous  pressures  are 
created  and  form  a  shock  front.  The  shock  front  propagates  approximately 
spherically  at  supersonic  speeds  and  produces  a  crusning  overpressure  force  with 
accompanying  gusts  of  dynamic  forces.  The  thermal  energy  effect  is  measured  in 
calories  per  square  centimeter  (cal/cm2),  and  blast  effects  in  pounds  per  square 
inch  (psi). 

Mathematical  models  are  available  to  study  the  thermal  and  blast  effects. 
Computer  codes  for  these  models  require  many  hours  of  computer  time;  consequently 
the  precise  codes  are  not  suited  for  systems  analysis  or  war  games.  Reliable 
models  based  on  the  precise  hydrodynamic  and  radiation  hydrodynamic  models  have 
been  developed,  tested,  and  improved  by  the  Air  Force  Weapons  Laboratory.  Codes 
for  these  models  require  only  milliseconds  of  computer  time  and,  therefore,  are 
suitable  for  systems  analysis.  The  computerized  versions  of  these  codes  bear 
the  names  SABER  and  SNAFT,  and  are  used  to  approximate  the  blast  and  thermal 
envi ronments ,  respecti vel y . 

The  Systems  Analysis  Blast  Environment  Routine  (SABER)  hcs  been  modified  to 
determine  only  the  ranges  of  given  levels  of  overpressure  and  the  times  of  shock 
arrival  at  those  ranges.  This  is  a  restricted  use  of  the  multipurpose  program. 
The  modified  version  is  called  SABERCM.  Inputs,  in  addition  to  specified  peak 
overpressure,  are  nuclear  weapon  yield,  height  of  burst,  terrain  height,  and 
aircraft  altitude.  Outputs  are  overpressure  range  and  time  of  shock  arrival. 

SNAPT  is  a  computerized  model  which  can  be  used  to  calculate  the  free-field 
thermal  energy  resulting  from  the  detonation  of  a  nuclear  weapon  or  to  calculate 
the  range  at  which  a  given  level  of  free-field  thermal  energy  occurs.  SNAPT  has 
been  modified  to  perform  only  the  latter  calculation  as  a  subroutine  named 
SNAPTCM.  Necessary  input  data,  other  than  the  free-field  energy  level,  consist 
of  nuclear  weapon  yield,  height  of  burst,  terrain  height,  aircraft  altitude,  and 
pertinent  atmospheric  conditions.  The  atmospheric  parameters  include  haze  layer 
height,  water  vapor  pressure,  ground  reflectance  (albedo),  and  visibility.  The 
horizontal  range  at  which  the  desired  free-field  thermal  energy  level  occurs  is 
output  from  SNAPTCM. 
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The  cookie-cutter  assumption  has  been  made  to  distinguish  between  regions 
of  lethality  and  nonlethality  to  aircraft.  For  these  purposes,  aircraft  vulner¬ 
ability  levels  for  the  thermal  and  overpressure  kill  mechanisms  are  specified 
in  cal/cm2  and  psi,  respectively.  Under  the  cookie-cutter  assumption,  aircraft 
are  assumed  to  be  killed  if  either  specified  vulnerability  level  is  exceeded,  . 
and  safe  otherwise. 

Vulnerability  levels  are  input  to  the  nuclear  routines  along  with  atmospheric 
conditions,  terrain  height,  and  weapon  characteristics  in  order  to  calculate  the 
lethal  nuclear  environmental  ranges  fof'  a  fixed  height  of  receiver.  A  point  of 
detonation  is  first  specified  for  each  missile  at  each  potential  target,  based 
upon  where  the  aircraft  from  that  airfield  are  located  when  the  missile  arrives. 
The  routines  SABERCM  and  SNAPTCM  require  a  height  of  receiver  to  compute  the 
horizontal  ranges  of  the  lethal  nuclear  environment.  The  height  of  receiver  is 
taken  as  the  altitude  of  the  aircraft,  accord’* ng  to  its  flight  profile,  at  the 
time  of  weapon  detonation  (relative  to  the  brake  release  time).  This  is  equiva¬ 
lent  to  slicing  the  spherically  propagating  shock  front  (or  overpressure  contour) 
and  the  thermal  contour  with  a  plane  parallel  to  the  ground  at  a  distance  above 
the  ground  equal  to  the  altitude  of  the  aircraft  at  the  above.  This  horizontal 
plane,  called  the  lethal  plane,  is  the  geometrical  structure  in  which  the  deter¬ 
mination  of  lethal  area  is  accomplished. 

The  general  appeararce  in  the  lethal  area  plane  of  the  overpressure  and 
thermal  contours,  at  the  lethal  levels  specified  by  the  vulnerability  levels  and 
relative  to  the  other  input  data,  is  that  of  two  concentric  circles  centered  at 
the  perpendicular  point  projection  of  the  burst  center  onto  the  plane.  Figure 
5  depicts  the  intersection  of  the  lethal  area  plane  with  the  letnal  overpressure 
contour.  The  horizontal  range  associated  with  the  lethal  overpressure  contour 
is  the  lethal  overpressure  radius,  and  is  measured  from  the  perpendicular  point 
projection  of  the  burst  center  onto  the  plane  to  the  lethal  overpressure  contour. 
Similarly,  the  lethal  thermal  radius  is  that  horizontal  range  associated  with 
the  lethal  thermal  contour.  Thus,  the  nuclear  routines  are  used  to  compute  the 
lethal  contours  needed  to  compute  lethal  area. 

The  orientation  of  the  aircraft  is  not  considered  in  computing  the  horizontal 
ranges  associated  with  the  overpressure  or  thermal  kill  mechanisms.  The  lethal 
thermal  radius  is  computed  under  the  assumption  that  the  aircraft  is  oriented 


20 


so  as  to  receive  the  maximum  amount  of  thermal  energy.  Similarly,  the  fact 
that  the  aircraft  is  better  equipped  to  withstand  tne  overpressure  shock  front 
in  one  position,  as  opposed  to  another,  is  also  not  taken  into  account  in  the 
computation  of  the  lethal  overpressure  radius. 

The  lethal  radii  which  are  outputs  from  the  nuclear  routines  are  computed 
for  a  stationary  receiver  (aircraft).  Thus,  the  lethal  contours  defined  by  the 
lethal  radii  must  be  adjusted  to  account  for  a  moving  aircraft.  This  is 
accomplished  in  the  lethal  area  plane. 

The  letha  'a  (within  the  lethal  area  plane)  resulting  from  the  detona¬ 
tion  of  a  nuc.  .  weapon  is  that  area  within  which  the  aircraft  cannot  survive 
if  located  there  at  the  onset  of  the  detonation.  Tabulated  data  for  the  air¬ 
craft  flight  profile  is  used  in  conjunction  with  the  distance  to  the  centroid 
(defined  in  section  V,  fifth  assumption)  to  transform  the  lethal  overpressure 
contour  in  the  lethal  area  plane  into  a  locus  of  points  describing  the  boundary 
of  the  area  reflecting  aircraft  kill  from  overpressure.  The  lethal  thermal 


21 


Ip 


k 


AFWL-TR- 73-242 


f 


i 


|  contour  remains  unchanged  since  the  thermal  energy  propagation  time  is  necli- 

|  gible.  The  entire  process  of  lethal  area  determination  takes  place  in  the 

|  lethal  area  plane. 
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A  top  view  of  a  typical  lethal  area  plane  containing  a  lethal  overpressure 
contour,  a  lethal  thermal  contour,  and  a  centroid  (designated  by  an  asterisk) 
is  offered  in  figure  6.  The  lethal  overpressure  contour  is  represented  by  the 
dashed-line  circle,  the  lethal  thermal  contour  by  the  solid-line  circle,  and 
the  perpendicular  point  projection  of  the  burst  center  onto  the  lethal  area 
plane  by  the  symbol  X.  The  remaining  task  is  to  adjust  the  overpressure  contour 
to  account  for  the  movement  of  the  aircraft. 
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Figure  6.  Top  View  of  Lethal  Area  Plane 
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Points  lying  on  the  lethal  overpressure  contour  are  used  to  determine  the 
boundary  for  the  overpressure  lethal  locus,  that  is,  the  lethal  area  at  the 
onset  of  the  burst  associated  with  the  overpressure  kill  mechanism.  The  air¬ 
craft  are  assumed  to  be  emanating  radially  from  a  point  called  the  centroid. 
The  points  on  the  lethal  overpressure  contour  are  backed  up  radially  toward 
the  centroid  by  the  distance  flown  between  detonation  and  the  arrival  of  the 
shock  wave  at  the  overpressure  contour.  This  distance  is  obtained  by  inter¬ 
polation  from  the  distance/time  coordinates  representing  the  aircraft  flight 
profile.  This  radial  translation  of  the  overpressure  contour  toward  the 
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centroid  usually  results  in  a  petal  or  egg-shaped  overpressure  lethal  locus  as 
shown  in  figure  7.  More  complex  shapes  may  result  when  the  centroid  is  within 
the  lethal  overpressure  contour. 


Figure  7.  Overpressure  Lethal  Locus 

The  overpressure  lethal  locus  encompasses  the  overpressure  kill  region  at 
the  time  of  detonation  under  the  cookie-cutter  assumption  for  a  moving  aircraft. 

An  underlying  assumption  is  that  the  aircraft  maintains  radial  flight  from  the 
centroid.  Since  the  aircraft  is  constrained  fay  the  aircraft  flight  profile,  if 
within  the  overpressure  lethal  locus  at  the  onset  of  the  burst,  it  will  be  inter¬ 
cepted  by  the  supersonically  propagating  shock  front  at  a  higher  level  of  over¬ 
pressure  than  it  can  withstand.  A  possibility  exists  that  the  aircraft  could  be 
located  within  the  lethal  thermal  contour,  as  well  as  within  the  overpressure 
lethal  locus,  at  the  onset  of  the  detonation. 

The  overpressure  lethal  locus  is  combined  with  the  lethal  thermal  contour  to  [ 
produce  the  boundary  of  the  lethal  area.  Figure  8  gives  an  example  of  a  lethal 
circle/petal  area  with  respect  to  a  moving  aircraft.  Numerical  integration  is  | 

I 

used  to  compute  the  area  within  this  lethal  region.  | 
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Figure  8.  Lethal  Area  for  Moving  Aircraft 

6.  DETERMINATION  OF  SURVIVABILITIES 

The  values  S..,  which  appear  in  the  objective  function  (maximizing  the 
expected  value  killed),  must  be  computed  for  each  weapon  in  group  j  against  each 
target  i.  Assuming  that  the  aircraft  at  each  target  are  uniformly  distributed 
over  some  area  at  each  weapon  arrival  time,  the  probability  of  kill,  Pk,  of 
each  aircraft  at  that  target  is  the  quotient  of  the  lethal  area  divided  by  the 
area  in  which  the  aircraft  could  be  located  when  the  weapon  arrives.  In  this 
report,  the  Pk  will  be  defined  as  the  probability  of  kill  given  that  the  weapon 
successfully  detonated  at  the  target.  In  QUANTO,  reliability  factors  are  given 
for  each  missile  type  for  the  probabilities  that  (1)  the  missile  is  successfully 
launched,  (2)  the  missile  successfully  reaches  the  target,  and  (3)  the  warhead 
successfully  detonates.  The  overall  reliability  of  a  missile  is  the  product  of 
these  three  reliabilities.  The  survivability  of  a  target  from  a  single  weapon 
is  than 

Survivability  =  1  -  P^  *  (reliability) 

The  P^  of  each  weapon  versus  each  target  is  computed  from  input  data  in 
QUANTO.  Initially,  the  arrival  time  of  the  missile  on  the  target  is  computed. 
Time  zero  is  the  time  at  which  all  of  the  first  missiles  from  each  submarine  are 
simultaneously  launched.  Subsequent  salvos  from  the  submarines  are  launched 
after  time  zero,  as  determined  by  the  salve  number  and  the  missile  launch 
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interval.  The  flight  time  of  a  missile  to  a  given  target  is  interpolated  from 
input  distance/time  missile  trajectory  information  after  the  distance  from  launch 
point  to  target  (i.e.,  the  coordinates  of  the  base)  is  computed.  The  arrival 
time  of  the  weapon  on  the  target  is  the  sum  of  the  launch  time  and  the  flight 
time. 


lhe  location  of  the  aircraft  at  the  time  of  weapon  arrival  may  be  computed  f 

from  the  input  aircraft  flight  profile  and  the  brake  release  times.  The  air-  ] 

craft  are  assumed  to  disperse  radially  from  a  single  point,  called  the  centroid.  | 

Unless  aircraft  are  assumed  to  be  departing  in  both  directions  from  a  base  (from  j 
dual  runways,  for  instance),  the  centroid  will  not  be  on  the  runway „  for  the  j  ; 

centroid's  location  is  a  function  of  the  time  it  takes  an  aircraft  to  raise  its  '  • 

gear  anc  flaps,  reach  a  turn  altitude,  etc.,  and  then  make  a  turn  to  its  fly-out  \  \ 

<  - 

direction.  The  distance  an  aircraft  will  be  from  the  centroid  at  weapon  arrival  ‘  : 
time  is  computed  by  (1)  subtracting  the  brake  release  time  from  the  weapon 
arrival  time  to  obtain  the  time  the  aircraft  has  had  to  escape  before  the  weapon 
arrives,  (2)  interpolating  in  the  aircraft  flight  profile  to  obtain  the  distance 
the  aircraft  has  traveled  from  brake  release,  and  (3)  subtracting  the  distance 
from  brake  release  to  centroid  from  the  total  distance  traveled.  The  distance 
to  t.ie  centroid  from  brake  release  point  is  input  for  each  target.  In  this 
manner,  QUANTO  computes: 


Rj  =  the  radial  distance  from  the  centroid  of  the 
first  aircraft  at  weapon  arrival  time 
and 

=  the  radial  distance  from  the  centroid  of  the 
last  aircraft  at  weapon  arrival  time 

If  an  aircraft  has  either  not  begun  its  takeoff  or  not  reached  the  centroid, 
its  radial  distance  from  the  centroid  is  set  to  zero.  The  intent  here  is  to 
treat  all  aircraft  which  have  not  reached  the  centroid  as  essentially  undis¬ 
persed  aircraft  which  can  be  targeted  with  a  single  SLBM. 

The  area  of  kill  generated  by  a  warhead  detonation  for  aircraft  of  a  given 
type  is  dependent  on  many  parameters,  as  described  in  the  preceding  section  of 
this  report.  Many  of  these  parameters  are  needed  to  describe  the  nuclear 
environment  and  are  directly  supplied  by  user  inputs.  First,  the  horizontal 
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ranges  are  determined  at  which  a  stationary  receiver  would  experience  a  lethal 
overpressure  or  thermal  effect.  Then  the  assumed  detonation  point  of  the  SLBM, 
the  aircraft  climb  profile,  and  the  distance  from  brake  release  point  to  the 
centroid  are  used  as  described  in  the  previous  section  to  determine  the  shape 
of  the  circle/petal  thermal /overpressure  lethal  area  at  detonation  time,  thus 
taking  into  account  the  moving  receiver  (aircraft).  The  lethal  area  varies 
somewhat  with  the  distance  of  the  detonation  from  the  centroid,  because  the 
aircraft  are  at  different  altitudes  and  velocities  at  different  points  in  the 
climb  profile.  Hence,  an  approximation  must  be  made  of  the  lethal  area  used  in 
the  calculation  of  P^,  and  QUANTO  must  make  some  assumption  about  where  the 
weapon  might  land  without  having  determined  yet  how  many  total  SLBMs  will  be 
allocated  to  the  target. 

QUANTO  assumes  that  the  attacker  can  compute  R  ,  assuming  a  certain  brake 
release  time,  and  accepts  this  as  the  farthest  distance  that  the  first  aircraft 
on  a  base  can  achieve  by  weapon  arrival  time.  However,  although  the  attacker 
might  also  be  able  to  compute  R^  based  on  the  stipulated  aircraft  takeoff 
intervals  and  brake  release  time,  he  would  realize  that  unanticipated  delays 
could  occur  (in  detection  and  warning  of  attack,  etc.)  and  might  consider  it 
equally  likely  to  find  aircraft  at  any  point  within  the  circle  of  radius  Rj. 
Consequently,  QUANTO  attacks  that  area,  in  general,  with  a  uniformly  dense 
distribution  of  weapons.  Thus,  an  average  weapon  might  land  at  a  distance 
(Rj//2T)  from  the  centroid  because  a  circle  of  radius  (R3//2T)  contains  half  the 
area  within  the  circle  of  radius  Rj.  One  benefit  of  this  uniform  attack  is 
that  a  del^v  in  the  brake  release  time  will  not  significantly  decrease  the 
expected  kill  and  may  result  in  a  large  increase  in  kill.  When  the  SLBM  deto¬ 
nates  at  distance  (R.//?)  from  the  centroid  and  the  entire  lethal  circle/petal 
area  falls  within  the  circle  of  radius  R1$  the  is  simply 

A^  [for  Rj//Srj 
k  nRx2 

where 

Al  [for  Rx//?j  (24) 

indicates  the  lethal  area  when  the  detonation  point  is  at  range  (R1/*//?)  from  the 
centroid.  This  situation  is  illustrated  in  figure  9. 
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Figure  9.  Detonation  Point  at  RjJ? 


The  use  of  equation  (24)  results  in  an  overestimate  of  the  expected  kill  of 
aircraft  in  the  case  of  very  few  aircraft  on  a  base  and  a  lethal  area  which  is 
a  large  portion  of  the  circular  area  nRj2.  With  few  (or  one)  aircraft,  R^  may 
be  only  slightly  less  than  (or  equal  to)  R}.  This  situation  is  shown  in  figure 
10  for  the  lethal  area  labeled  A.  In  this  case,  equation  (24)  predicts  a  large 
percentage  of  the  aircraft  killed,  although  when  the  brake  release  time  is 
certain,  no  kills  result.  To  guard  against  this  possibility,  the  attacker 
would  wish  to  reduce  his  estimate  of  Pk  in  allocating  his  attacking  weapons. 

One  way  in  which  he  could  do  this  would  be  to  replace  the  procedure  of  the 
previous  paragraph  with  one  in  which  the  weapon  was  placed  at  distance 
(Ri  +  Rn)/2  from  the  centroid,  assuming  the  aircraft  were  uniformly  distributed 
throughout  the  annulus  of  radii  Rj  and  R^.  This  Pk  is  the  shaded  area  within 
the  circle/petal  labeled  B  in  figure  10,  divided  by  the  annulus  area  ji^Rj2  -  RN2j 
The  shaded  area  within  B  may  be  approximated  by  considering  the  circle/petal  B 
as  a  circle  of  equivalent  area  centered  at  (Rj  *  Rj^)/2  radial  distance  from  the 
centroid,  and  computing  the  area  within  both  the  equivalent  circle  and  the 
annulus  thickness.  This  latter  common  area,  which  will  be  labeled  jfor 
(Rj  +  RN)/2]  ,  may  be  computed  from  closed-form  geometric  expressions.  The 
single  aircraft  situation  is  handled  in  this  manner  by  artificially  setting 
RN  =  Rx  -  0.01  so  that  the  annulus  has  a  small  positive  area.  The  Pk  formula 
thus  becomes 


Pk  =  min 


Al  [for  R,//?]  Aun  [for  (R2  +  RN)/2* 


Jl  Rj2 


n  R,»  - 


(24) 
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Figure  10.  Lethal  Area  Detonation  Illustration 


The  use  of  this  formula,  when  the  second  quotient  is  the  minimum  of  the  two, 
encourages  OUANTO  to  allocate  a  second  weapon  to  the  target  since  the  first 
weapon's  is  lowered.* 

As  shown  in  the  previous  figures,  the  thermal  lethal  circle  usually  extends 
farther  from  the  centroid  than  does  the  overpressure  lethal  petal.  Likewise, 
when  the  SLSM  detonates  at  the  centroid,  the  thermal  circle  usually  extends 
farther  from  the  centroid  than  does  the  region  of  overpressure  kill,  which  in 
this  special  case  is  also  a  circle.  "Usually"  in  this  context  means  for  those 
combinations  of  flight  profile  and  overpressure/thermal  vulnerability  levels 
normally  of  interest.  However,  for  some  hardness  levels  and  aircraft  flight 
profiles,  the  region  of  overpressure  kill  may  totally  encompass  the  thermal 
lethal  circle. 

The  farthest,  reach  of  the  lethal  region  in  a  direction  away  from  the  centroid 
becomes  a  concern  when  is  small  enough  that  the  lethal  area  may  protrude 
beyond  R.  for  a  given  weapon  placement.  When  the  lethal  area  so  protrudes,  the 


*The  validity  of  the  met.iod  implicit  in  the  second  term,  and  indeed  of  the 
whole  procedure  in  equation  (24),  has  been  confirmed  by  comparing  answers 
obtained  from  QUANTO  with  those  obtained  using  a  Monte  Carlo  simulation 
model . 
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first  expression  competing  for  the  minimum  in  equation  (24)  is  in  error  because 
the  aircraft  cannot  be  located  in  the  protruding  portion  of  the  lethal  area. 

To  handle  these  cases,  the  weapon  is  assumed  to  detonate  at  a  position  for 
which  no  protrusion  occurs,  if  possible. 

Therefore,  QUANTO  computes 

Rlr  [for  R1//2']=  the  distance  from  detonation  point  to  lethal 
region  boundary,  in  a  direction  away  from  the 
centroid,  for  a  detonation  at  distance 
(R^/?)  from  the  centroid. 

It  will  usually  be  true  that 

RLR  [for  Ri  "  RLR  Lfor  “  RLR  [for 

i.e.,  R^r  varies  little  as  the  detonation  point  ~  adjusted  to  avoid  protrusion. 
Now,  if 

R1  <.  R]_R  [for  0],  then  QUANTO  sets  =  1 


But  if 


Ri  >  RLR  *-for 


(R1//^)  +  Rlr  [for  R1//2’]>  Rj  (i.e.,  the  lethal  area 

protrudes  beyond  Rx) 

thtn  QUANTO  computes  the  as 

Al  [for  Rj  -  RLR  [for  Rj/^]J  ALAN  [for  (Rx  +  RN)/2] 

p  =  min  - ,  - - - - - 

11  Ri2  n(Ri2-RN2) 


n  R:2 


where  the  first  quotient  assumes  placement  of  the  weapon  at 


(*i  -  \r  [for  V^J) 


in  order  to  avoid  protrusion. 
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Thus,  the  complete  formula  for  is 
if  R1  ^  Rlr  [for  0] 

Al  [for  Rj  -  Rlr  [for  R,//?]!  AUN  [for  (R,  +  RN)/2] 

"  Ri2  ’  n(R,2  -  RN2) 

if  Rx  >  R^r  [for  0]  and  Rj/v^T  +  R^R  [for  Rx/  T]  >  Rx 

Al  [for  Rj//?]  Alan  Jfor  (Rx  +  Rn)/2J  J 

- ,  - - - - - \  otherwise  (26) 

"  R,2  n(Ri2  -  rn2)  ) 

It  may  be  noted  in  the  above  formula  that  the  weapon  placement  may  be  at  dis¬ 
tances  from  the  centroid  of  0  (when  P^  =  1),  Rx  -  R^R  [for  R^*^],  (Rx  +  RN)/2, 
or  R^/?.  This  will  be  of  concern  when  more  than  one  aircraft  type  is  considered 
in  the  model. 

It  should  be  noted  that  these  formulas  for  Pk  are  inaccurate  when  the  dis¬ 
tribution  of  aircraft  is  far  from  uniform  over  an  area,  as  might  be  the  case  if 
a  number  of  aircraft  had  not  left  the  base  by  the  time  of  a  weapon  arrival. 

This  situation  will  not  occur  if  the  aircraft  beddown  is  a  rational  one,  intended 
to  prevent  mass  kills  by  single  weapons.  Because  of  the  assumption  of  uniform 
distribution  of  aircraft,  QUANTO  will  underestimate  the  aircraft  kills  in  these 
situations.  However,  a  simulation  program  may  be  used  to  discover  if  such 
conditions  exist  and  to  estimate  the  resultant  kills. 

If  the  aircraft  were  actually  uniformly  distributed  over  the  areas  assumed, 
the  actual  P^  values  (and,  hence,  S^j  values)  realized  by  the  allocated  weapons 
would  not  agree  exactly  with  those  computed  by  formulas  (23)  and  (24).  This  is 
due  to  the  impossibility  of  determining  the  realized  lethal  area  sizes  before 
determining  the  number  of  weapons  (and,  therefore,  the  precise  placement  of 
weapons)  on  each  target.  The  computed  weapon  allocation  is  optimal  for  the  P^ 
values  computed.  However,  the  plus  and  minus  errors  between  realized  and 
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computed  P^'s  for  individual  weapons  tend  to  balance  out  to  a  small  overall 
error  when  summed  over  all  the  weapons.  This  is  because  the  computations  of 
Pk's  are  based  on  an  average  placement  or  each  weapon  on  each  target. 

It  has  been  assumed  that  the  survivabilities  S. ^  are  independent.  Thus,  no 

*  J 

weapon  on  target  i  can  cause  collateral  damage  on  aircraft  from  another  base, 
and  the  area  purged  of  aircraft  by  a  detonation  can  become  populated  to  an  equal 
aircraft  density  by  other  aircraft  before  the  next  weapon  arrives.  Detailed 
base-by-base  simulation  of  the  attacks  produced  by  QUANTO  has  shown  that  the 
actual  resultant  kill  (output  from  the  simulator)  is  not  significantly  different 
from  QUANTO's  predicted  kill. 
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SECTION  IV 

THE  COMPUTER  PROGRAM 

Figure  11  shows,  in  general  terms,  the  operation  of  QUANTO.  After  the  data 
for  the  problem  is  input,  areas  of  lethality  and  the  survivabilities  $. .  are 

*  J 

computed.  The  optimal  missile  laydown  n.  -  is  then  determined.  QUANTO  will  then 

'  J 

relocate  a  submarine  to  a  better  position,  if  the  user  has  requested  submarine 
optimization,  and  recompute  the  optimal  n- -  for  the  new  positions  of  the  sub- 

*  J 

marines,  'fter  submarine  optimization  is  completed,  aircraft  may  be  relocated 
to  improve  the  number  surviving,  with  n-  •  recomputed  following  each  shift  of 

*  J 

aircraft.  When  beddown  has  been  optimized,  the  optimal  nonintegral  V.  and  n.j 
are  integerized  and  final  output  is  produced. 

Figure  12  indicates  several  additional  details  of  QUANTO.  The  criteria  for 
terminating  submarine  optimization  and  beddown  optimization  are  indicated  in 
test  blocks.  A  mode  parameter,  input  on  the  first  data  card  of  a  problem  deck, 
controls  from  whence  input  is  taken  and  how  much  of  the  program  is  executed. 
Table  II  described  the  mode  options.  These  options  permit  the  user  to  observe 
partial  computations  tor  validity  without  risking  a  large  expenditure  of 
computer  time. 

The  principal  subroutines  .  •  QUANTO  (QUANTO  being  the  name  of  the  main 
program)  and  their  functions  are  listed  in  table  III. 


OUTPUT  EXPECTED 
KILL 


COMPUTE  INITIAL  KILL  AND  MULTIPLIERS 


COMPUTE  KILL  AND  OUTPUT 


*■ m 


Table  II 
MODE  OPTIONS 

I 

Mode  QUANTO  Operations 

0  Input  all  data  except  convergence  parameters  and  initial 

allocation  n. .  from  cards. 

•  J 

Terminate  problem  after  computing  nuclear  effects  and 
constructing  flight  profiles. 

1  Input  all  data  from  cards. 

Terminate  after  computing  survivabilities  S..  and  kills 

*  J 

resulting  from  initial  allocation. 

Write  information  for  automatic  program  restart  on  tape. 

2  Input  all  data  (except  parameters  on  first  card)  from 

restart  tape. 

Terminate  problem  computations  prior  to  completion  only 
if  the  time  limit  for  processing  is  reached,  at 
which  time  a  restart  tape  will  be  written. 

3  Input  all  data  from  cards. 

Terminate  problem  as  for  mode  =  2. 


* 
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Routine 

Name 


QUANTO 


Table  III 

PRINCIPAL  ROUTINES  AND  THEIR  FUNCTIONS 

Description 

Main  Program. 

QUANTO  reads  the  input  data,  either  from  cards  or  from  a  restart 
tape,  as  controlled  by  the  mode  parameter  on  the  first  input 
card  of  each  problem  deck.  The  data  describing  the  problem 
are  printed,  and  if  the  data  are  read  from  cards,  a  summary  of 
the  input  data  is  also  printed.  The  aircraft  profile  and 
parameters  affecting  nuclear  effects  are  not  read  by  QUANTO, 
but  by  PROCESS,  called  by  QUANTO. 

Computations  of  survivabilities  S.  .  are  performed  mainly  ;n  QUANTO. 

•  J 

DETAREA  provides  QUANTO  the  necessary  lethal  areas,  but  QUANTO 
computes  flight  times  and  distances  and  the  resultant  set  of 
Pt  and  S. .  values,  with  the  help  of  interpolation,  look-up, 
and  distance  computational  routines. 

After  input  of  the  initial  allocation,  QUANTO  controls  the  sequenc¬ 
ing  of  operations  in  the  iterative  procedure  for  optimizing  the 
missile  laydown.  When  provided  an  by  ADJLAM,  OUANTO  recomputes 
n..,  S  -  - ,  x..,  and  f(n,-i).  Control  of  the  iteration  cutoff  and 

1 J  I J  *  J  1 

intermediate  output  is  accomplished  in  QUANTO. 

Relocation  of  aircraft  is  completely  performed  in  QUANTO,  but 
QUANTO  calls  other  routines  for  submarine  optimization  and 
Integerization  of  laydown  and  beddown. 
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Table  III  (cont'd) 

Routine 

Name  Description 

PROCESS  Called  by  QUANTO. 

PROCESS  reads  aircraft  profile  data  and  nuclear  effects  parameters 
from  cards  for  each  type  of  aircraft.  PROCESS  generates  dis¬ 
tance/time  coordinates  for  each  aircraft  for  the  specific 
altitude  of  level-off.  The  input  data,  as  well  as  the  generated 
distance/time  coordinates,  are  output.  If  the  lethal  over¬ 
pressure  and/or  thermal  radius  and  the  time  of  shock  arrival  are 
not  present  in  the  input,  PROCESS  computes  these  values  for  the 
yield  of  each  type  of  missile.  The  nuclear  effects  information 
is  summarized  in  the  output  from  PROCESS. 

DETAREA  Called  by  QUANTO. 

From  given  aircraft  profiles  and  geometry  of  flyout  and  detonation, 
DETAREA  computes  the  lethal  area  with  respect  to  the  moving  air¬ 
craft,  i.e.,  the  circle/petal  area  describing  the  thermal/ 
overpressure  kill  region. 

ALOUT  Called  by  QUANTO. 

ALOUT  produces  a  list  of  the  allocation  n.-  in  two  formats.  First, 

*  J 

by  target:  the  missiles  allocated  to  that  target  are  listed  in 
order  by  submarine  number  and  salvo  number  within  the  submarine. 
Second,  by  submarine:  the  missiles  are  listed  in  order  by 
salvo,  together  with  the  targets  to  which  they  are  allocated. 
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Table  III  (cont'd) 


Routine 

Name  Description 


ALIN T  Called  by  QUANTO. 

AUNT  integerizes  the  allocation  matrix  n-  -.  This  process  is  not  a 

simple  rounding  of  the  nonintegral  n • .  values*  but  an  assignment 

*  0 

of  integral  values  to  the  highest  fractional  parts  so  as  to 
M 

make  z  n- -  =  N.  for  each  j. 
i=j  J 

VINT  Called  by  QUANTO. 

VINT  integerizes  the  beddown  values  V-  in  a  manner  similar  to  AUNT. 
TGTKIL  Called  by  QUANTO. 

L  "ii  , 

TGTKIL  computes  the  survivability  products  n  S- •  J  (for  each 

0=1 

aircraft  type),  the  x,.  values,  the  number  of  aircraft  killed 

'  J 

at  each  base,  and  a  rough  idea  (obtained  by  rounding  n-.  values) 

*  J 

of  the  number  of  weapons  allocated  to  each  base. 


ADJLAM  Called  by  QUANTO. 


ADJLAM  finds  An  by  first  finding 

Akfc  =  m1n  lAig,  suc^  ni£  a  0*0001|  i  =  1,  2,  ...,  M,  and 

V  =  “*  !xnf  1  “  ’• 2 . "• such  that  lk( <  V  - E 

for  some  weapon  group  i  and  some  tolerance  e.  Then  ADJLAM 
either  computes  the  proper  allocation  adjustment  An  to  force 
the  value  of  xk£  toward  the  value  of  xm£  (for  the  single  t>pe 
aircraft  model)  or  calls  the  function  XNEWT  to  compute  An  (for 
the  mixed  force  model).  Tallies  of  weapons  allocated  to  each 
base  are  updated  after  the  change  of  n-.  by  the  An  adjustments. 

•  J 
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Table  III  (cont'd) 

Routine 

Name  Description 

SUBADJ  Called  by  QUANTO. 

SUBADJ  locates  the  least  effective  submarine  and  the  submarine 
location  with  the  most  potential,  as  described  in  the  test  on 
submarine  optimization.  The  submarine  is  relocated  to  t  e 
better  position  and  its  missiles  are  allocated  to  bases  having 
high  X . .  values,  as  described  in  the  text. 

XAREA  Called  by  QUANTO. 

XAREA  computes  the  area  of  intersection  of  a  circle  and  an  annulus, 
under  all  conditions  of  annulus  radii,  circle  radius,  and 
offset  of  circle  center. 
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SECTION  V 
ASSUMPTIONS 

In  this  section,  :he  principal  assumptions  in  QUANTO  are  described  and  are 
briefly  discussed. 

Assumption:  All  input  SLBMs  are  used  against  aircraft,  i.e.,  the  attacking 
force  decides  what  portion  of  its  SLBMs  to  use  against  the  flushing  aircraft 
force  prior  to  running  a  problem  and  the  SLBMs  in  QUANTO  represent  that  portion. 
Of  course,  for  each  submarine,  only  a  partial  load  of  missiles  need  be  input  for 
a  problem. 

Assumption:  The  survival  probabilities,  S. ...  are  independent.  No  collateral 

*  J 

damage  may  affect  an  aircraft  departing  one  base  as  a  result  of  a  detonation  of 
an  SLBM  allocated  to  another  base.  Furthermore,  the  effectiveness  of  a  weapon 
in  group  j  on  base  i  is  measured  by  $..,  a  numerical  quantity  which  is  independ- 

'  J 

ent  of  the  number  of  weapons  which  have  previously  arrived  or  will  subsequently 
arrive.  Stated  differently,  the  area  in  which  aircraft  may  be  located  at  the 

time  of  a  later  weapon  arrival  is  not  considered  to  contain  voids  left  by  pre¬ 

viously  arriving  weapons.  Detailed  Monte  Carlo  simulation  of  QUANTO-produced 
attacks  shows  that  QUANTO's  predicted  kill  (using  such  S-.  values)  is  close  to 

*  J 

the  actual  kill  resulting  from  the  simulation. 

Assumption:  Thermal  and  overpressure  effects  have  lethality  according  to  a 
cookie-cutter  criterion.  In  other  words,  an  ai-craft  with  hardness  of  x  psi  and 
y  cal/cm2  is  killed  if  it  experiences  either  of  these  levels  or  higher,  but  is 

safe  from  (x-e)  psi  and  (y-e)  cal/cm2  for  any  e  >  0,  no  matter  how  small. 

Assumpti on :  At  all  times,  aircraft  are  uniformly  distributed  within  a  maxi¬ 
mum  circte,  defined  by  the  first  aircraft's  range,  the  area  of  which  is  contin¬ 
ually  increasing  with  time.  Thus,  the  survivabilities  S.j  are  computed  assuming 
the  attacker  will  pattern  his  weapons  for  uniform  coverage  of  the  maximum  circle 
of  aircraft.  In  some  cases  (few  aircraft  at  early  weapon  arrival  times),  this 
assumption  is  modified  to  allow  computation  of  S^.  by  assuming  instead  that  the 
aircraft  are  uniformly  distributed  throughout  an  annulus.  In  this  way,  a  lower 
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computed  S. .  results  and  more  realistic  expected  kills  result.  Of  course,  if 

*  *J 

the  aircraft  do  not  disperse  in  a  circular  pattern,  and  the  attacker  were  to  be 
granted  advance  knowledge  of  these  flyout  tactics,  a  greater  expected  kill  would 
result  since  the  weapons  have  a  smaller  area  to  attack. 

Assumption:  The  aircraft  radially  emanate  from  a  point  called  the  centroid 
of  the  aircraft  which,  for  a  given  flyout  profile  and  turn  geometry,  may  be 

determined.  The  distance  of  the  ce  roid  from  brake  release  point  is  a  param-  ! 

eter  which  may  be  input  for  each  ba--.~ ,  it  is  based  on  the  distance  the  aircraft  ; 

flies  without  turning  while  raising  its  gear,  climbing  to  turn  altitude,  etc. 

Assumption:  In  the  computation  of  S- .,  the  detonation  point  of  the  weapon  { 

1  J  5 

is  at  one  of  several  places,  as  described  in  another  portion  of  this  report.  j 

The  assumption  of  detonation  point  is  such  as  to  be  in  agreement  with  the  j 

uniform-attack-of-the-aircraft-area  assumption,  with  a  modification  of  the  j 

location  (1)  when  protrusion  of  the  lethal  area  beyond  the  maximum  circle  occurs,  j 

(2)  when  an  annular  P.  computation  yields  a  better  estimate  of  S..,  or  (3)  when  1 

k  1  j  \ 

a  weapon  on  the  centroid  kills  all  aircraft  of  a  single  type.  In  this  way,  some  *  1 

pains  are  taken  to  compute  S..  based  on  a  reasonable  estimate  of  the  weapon  i 

•  J  ; 

location,  without  knowledge  of  where  other  weapons  are  allocated.  .  ■ 

Assumption:  When  multiple  aircraft  types  are  included  in  the  model,  all  1 

aircraft  radially  emanate  from  a  single  centroid  and  weapons  are  patterned  to 
attack  uniformly  the  area  of  all  aircraft  types.  Point  values  of  aircraft  of 
different  types  may  make  some  bases  more  attractive  than  others;  but  on  those 
bases,  the  attack  is  assumed  to  be  uniform.  i 


j 
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SECTION  VI 

LAGRANGE  MULTIPLIERS  IN  THE  MIXED  FORCE  ALLOCATION  PROBLEM 


When  more  than  one  type  of  aricraft  may  be  leaving  each  base,  the  problem 
of  determining  the  optimal  missile  laydown  is  more  complicated  than  the  pre¬ 
viously  described  model.  For  purposes  of  explanation,  the  following  will  assume 
two  aircraft  types,  bombers  and  tankers,  indicated  by  B  and  T  subscripts, 
respectively.  However,  the  procedures  are  general  for  any  number  of  aircraft 
types.  The  objective  to  be  maximized  in  the  mixed  force  allocation  problem  is 


M  f  L  n.  " 

f(nij)  =  E  viB  i  -  s.jB  3 


with  the  same  stockpile  constraints 


M 

+  £  v< 


iT 


L  n. . 

i  -  n  s..T  13 

.  j=i  131  . 


(27) 


M 

Ti  n,*-i  =  N.,  i  =  1  *  2,  . . . ,  L 
j=l  13  3 

The  function  f(n-jj)  is  now  the  expected  kill  value  of  both  bombers  and  tankers 
with  a  sum  over  the  target  index  i  for  each  type  of  aircraft.  The  values  of 
the  bombers  and  tankers,  respectively,  leaving  base  i  are  V^B  and  V.y,  where 
each  type  aircraft  may  be  worth  a  different  amount  of  value  per  aircraft.  The 
survivabilities  S.^.g  and  S.^j  of  the  bombers  and  tankers,  respectively,  of 
target  i  from  a  weapon  in  group  j,  must  be  computed  slightly  differently  than 
the  previous  S... 

'  J 

Since  the  types  of  aircraft  may  have  different  thermal  and  overpressure 
hardness  levels,  there  is  a  circle/petal  combination  for  each  aircraft  type  at 
a  single  weapon  detonation  point.  Using  the  same  notations  as  before,  with  the 
additional  B  (bomber)  and  T  (tankers)  subscripts  (following  the  slashes),  the 
following  formulas  for  Pk^B  and  P^y  are  used  for  the  computation  of  bomber  and 
tanker  P^'s. 

In  the  mixed  force  allocation  problem,  it  will  be  assumed  that  the  SLBMs  are 
aimed  uniformly  at  the  entire  area  of  all  aircraft  (of  all  types).  When  the 
geometry  is  such  that  the  tankers  are  within  the  bombers,  as  shown  in  figure  13, 
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Figure  13.  Bomber/Tanker  Mix 


the  P^T  may  be  approximated  by  the  product  of  (1)  the  probability  that  a  random 
placement  wUhin  the  bomber  circle  lands  within  the  tanker  circle,  and  (2)  the 
probability  of  kill  of  a  tanker  given  that  the  weapon  detonates  within  the 
tanker  circle.  If  A^y  represents  the  lethal  area  of  the  SLBM  against  tankers, 
the  product  is 


P  -  nRT  x  AL/T  =  \/T 
k^T  nR|  jir|  nR| 


where  overlap  of  the  lethal  area  over  the  circle  of  radius  Ry  has  been  ignored. 
Thus,  for  these  assumptions  and  approximations,  P^y  is  independent  of  Ry. 
Should  the  bombers  be  within  the  tankers. 


_  AL/B 


so  the  maximum  radius  always  appears  in  the  denominator. 
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For  the  mixed  aircraft  model,  let  R2  and  R^  represent  the  radial  distances 
of  the  most  and  least  distant  aircraft,  respectively,  from  the  single  centroid 
(the  same  point  for  all  aircraft  types)  at  the  time  of  weapon  arrival.  The 
subscript  MAX  means  the  largest  for  all  aircraft  types,  and  other  notations  are 
analogous  to  those  used  in  the  discussion  of  the  single  type  aircraft  model. 

A.  If 


Ri  -  RLR/MAX  tfor 

then 


and 


assuming  that  the  weapon  detonates  at  the  centroid. 


B.  If 


(30) 


(31) 


(32) 


Ri  >  RLR/MAX  ffor 
and 


Rj/*^  +  [for  Rj/vT]  s.  Rj 


(33) 


so  that  neither  circle/petal  lethal  area  protrudes  beyond  R  ,  then 

[^or  Ri/^] 


0)  Pk/T  =  min 


1, 


JI  R, 


(34) 


and 


Pk/B  -  min 


\/B  Cfor  M'5’] 


n  R, 


(35) 
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if  either 


or 


Otherwi se , 


and 


alan/t 

for 

(R.  +  rnH 

\/T  IV 

ni 

(Ri2 

•  rn2) 

n  Rj2 

alan/b  I 

for 

(R1 +  %)'2] 

AL/B  tfor 

n! 

(R.2 

-rn2) 

n  R:2 

/o  \  r* 

alan/t  [for  (R1  +  rn)/2] 

(2)  P 


k/T 


k/B 


n(Ri2  “  rn2) 

alan/b  [f0r  (Ri  +  rn)/2] 

n(Ri2  -  RN2) 


(36) 


(37) 


(38) 


(39) 


In  words,  if  the  annulus  computations  in  equation  (38)  above  are  both 
smaller  than  their  corresponding  circular  Pk  computations  in  equation  (34), 
then  the  annulus  P^  formulas  are  used  for  all  aircraft  types.  Otherwise,  the 
formulas  in  equation  (34)  are  used. 

C.  Finally,  if  either  circle/petal  protrudes  when  positioned  at  Rj//?, 
the  weapon  is  moved  toward  the  centroid  as  before  and 


ind 


AL/T  [tcr  Ri  '  RLR/KAX  ^for  V^] 


AL/B  [for  Ri  "  RLR/MAX  ^for 

?k/B  =  TrT 


if  either 


(40) 


(41) 
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ALAN/T  [for  (Ri  +  rn)/2]  AL/T  [for  Rl  '  RLR/MAX  Lfor 


n  Rx2 


or 


"(R12  -  RN2) 

alan/b  [for  (R;  +  rn)/2]  al/t  [for  Ri  '  rlr/max  Cfor 


(42) 


"(Ri2  '  rn2) 


Otherwise, 


(2)  P 


and 


n  Rj2 


'YAN/T  [for  (Ri  +  Rn)/2] 

k/T  ?(*i2  -  rn2) 

alan/b  [for  (R1  +  rn)/2] 


(43) 


(44) 


IVB 


"(R.2  -  RN2) 


(45) 


It  should  be  noted  that  these  formulas  have  been  written  in  a  slightly  different 
form  to  ensure  that  the  Pk  computations  for  different  aircraft  types  are  all 
based  on  the  same  placement  of  the  weapon. 

The  survivabilities  are  simply 


SijB  =  1  "  Pk/B  *  (re11ab11ity  of  weaPt>n  ln  group  j) 


-  P^y  *  (reliability  of  weapon  in  group  j) 


The  technique  for  solution  of  the  constrained  maximization  problem  for  the 
mixed  aircraft  force  is  very  similar  to  the  techniques  for  the  previous  model. 
The  Lagrangian  function  for  the  new  objective  function  f(njj)  is 


(46) 


Setting  the  partial  derivatives  equal  to  zero,  as  before,  yields 


J* L«  x  +  x 

3nk£  'kfcB  Ak*T 


+  x. 


=  0 


k  =  1,  2, 


M;  l  =  1,  2,  ...»  L 


(47) 
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where 

Xk£B  ~  "VkB  {zn  Sk£B)  A  Skj3  ^ 

v  * 

and 

Xk£T  "  "VkT  (£n  Sk£l)  A  SkjT  kj 

J  ' 

For  the  mixed  aircraft  model,  the  new  definition  of  is 

Xk£  5  xk£B  *  Xk£T 


(48) 


(49) 


(50) 


or  in  the  case  of  more  than  two  aircraft  types,  X^  would  be  the  sum  of  the 
lambdas  corresponding  to  each  aircraft  type.  Fixing  £  and  letting  k  vary 
results  in  the  system  of  equations 


Xk£  ~  _X£*  k  ~  2.  ...»  M 

which  has  the  same  appearance  as  in  the  single  aircraft  model,  although  xk?  is 
differently  defined. 

The  iterative  procedure  is  again  based  on  finding  Xk£  <  Xm£  with  n^  >.  0.0001 , 
and  choosing  An  so  that  the  new  values  of  x^^  and  x^,  say,  xk?  and  >  ,  become 
equal.  This  value  of  An  is  the  root  of 


g(An)  =  ;;k£  -  x 


m£ 


"  xk£B  +  xk£T  '  xm£B  "  xm£T 

-  c  ~An  .  .  c  -An  .  c  +An  .  p  +An  .  /ci  \ 

'  5:<£B  Ak£B  ^k£T  Ak£T  m£B  Am£B  TotT  Am£l  1 

The  root  of  An*  of  g(Ao)  is  found  by  Newton's  successive  approximation  method, 
where 


gUn-j) 

Ani+i  =  a^.  -  gT^y  ,i=l,2 . 


(52) 
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In  QUANTO,  when  lAn^-j  -  An.  |  <  e  (the  same  tolerance  used  in  the  test  for 
convergence  of  and  *m£).  An*  is  set  equal  to  Ani+1.  This  iterative  formula 
may  diverge  (as  tested  by  |An.j+.|i  >.  200  in  QUANTO)  for  a  given  selection  of  AnQ. 
Convergence  is  attempted  for  AnQ  selections  of  n^,  0,  and  nk£/2,  successively, 
until  the  iteration  successively  converges.  If  the  root  An*  >  n^.  An  is  chosen 


as  n.,,;  i.e.. 


An  =  min  |An*,  nkJJ 


so  as  to  keep  (nk£  -  An)  nonnegative. 

This  choice  of  An  results  in  the  maximum  increase  in  the  expected  kill.  The 
proof  of  this  statement  is  quite  similar  to  the  analogous  proof  in  the  single 
aircraft  case.  If  the  kill  contribution  to  the  objective  function  f(n-jj)  from 
targets  k  and  m  after  An  weapons  are  moved  from  target  k  to  target  m  is  indicated 
as  r(An),  then  the  unconstrained  maximum  of  r(An)  occurs  where  r'(An)  =  0.  But 
this  is  at  An*,  the  root  of  g(An),  since 

r'(An)  =  -g(An)  (54) 

The  maximum  of  r(An),  constrained  by  0  <.  An  <  nt  ,  occurs  at  An  =  nt  if  An*  > 


nk£,  since 


I"<°)  ■  Sn,*  -  *k,  *  0 


r"  «  -s' (in)  <  0 


The  ‘terative  procedure  for  the  mixed  aircraft  model  is  exactly  the  same  as  that 
for  the  single  aircraft  model,  with  the  exception  of  the  new  definition  of  Akjl 
and  the  Newton  procedure  for  finding  An*.  The  submarine  relocation  procedure  is 
also  based  on  the  new  =  x^g  +  xk£j. 

The  beddown  optimization  procedure  is  slightly  changed  in  that  each  beddown 
change  simultaneously  moves  some  of  each  type  of  aircraft.  Thus,  for  bombers, 
the  value 


6Vb  ■  v(j»,  VjB%a  -  smki) 
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is  shifted  from  base  k  to  base  m  where  these  bases  have,  respectively,  the 
lowest  and  highest  survivability  products 


n  s..R 

j=i  ljB 


ij 


but  only  those  bases  k  for  which  V^g  >  0.0001  compete  for  the  lowest  product. 
The  value  of  tankers  shifted,  aVj,  (computed  like  L\Ut  replacing  "B"  subscripts 
with  "T"  subscripts)  depends  on  the  tanker  survivabilities  S^j  and  values  V-T, 
and  thus,  bomber  and  tanker  relocations  may  involve  different  pairs  of  bases. 
The  beddown  optimization  stops  when  the  total  number  of  aircraft  (of  all  types) 
to  be  moved  in  a  single  beddown  change  is  less  than  0.05. 
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SECTION  VII 

BENEFITS  AND  FUTURE  USES  OF  QUANTO 

QUANTO  was  developed  to  investigate  the  sensitivity  of  total  bomber  force 
survivability  to  variations  in  thermal  and  overpressure  hardness.  In  the  process 
of  this  development,  numerous  other  parameters  have  been  included  as  variables 
in  the  model.  Consequently,  the  model  is  useful  for  evaluating  the  sensitivity 
of  surviving  aircraft  to  changes  in  aircraft  beddown  and  flight  profiles,  numbers 
and  types  of  submarines  or  missiles,  SLBM  performance  characteristics,  and  reac¬ 
tion  times,  as  well  as  aircraft  hardness.  Alternative  missile  laydowns,  sub¬ 
marine  locations,  and  aircraft  beddowns  may  be  compared  and  evaluated  using 
QUANTO  to  compute  expected  kills.  Contractual  studies  may  be  evaluated  for 
validity  and  contrasted  with  QUANTO  to  aid  in  understanding  their  results.  In- 
house  and  intra-AF  investigations  are  facilitated  by  the  availability  of  QUANTO. 

QUANTO  has  a  great  deal  of  flexibility.  It  is  relatively  fast  and  easy  to 
use  compared  to  other  flush  models.  An  optimal  laydown  may  be  computed  in  1  to 
3  minutes  of  computer  time,  submarines  may  be  optimized  in  about  5  minutes,  and 
optimal  beddowns  require  up  to  10  minutes,  where  these  times  are  largely  depend¬ 
ent  on  the  quality  of  the  selection  of  initial  laydowns,  submarine  positions, 
and  beddowns.  Each  weapon  is  considered  a  separate  entity,  not  as  a  member  of 
one  of  a  fixed  set  of  predefined  patterns.  In  addition,  QUANTO  permits  multiple 
types  of  aircraft  and  SLBMs,  each  with  its  own  performance  characteri stirs.  The 
modular  construction  of  QUANTO  permits  investigation  of  selective  changes  in  the 
assumptions  upon  which  the  model  is  based,  with  selective  program  changes. 

The  projected  future  uses  of  QUANTO  include  the  evaluation  of  other  models 
and  results  of  flush  studies,  studies  of  the  effects  of  parametric  variations 
on  the  survivability  of  a  mixed  force,  and  in-house  experimentation  and 
sensitivity  analysis. 
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APPENDIX  I 

APPLICATIONS  OF  LAGRANGE  MULTIPLIERS 
The  Basic  Problem 

The  computer  program  called  QUANTO  uses  the  Lagrange  multiplier  method  to 
optimize  the  allocation  of  weapons.  This  appendix  is  provided  to  introduce  the 
reader  to  the  basic  fundamentals  of  the  technique. 

The  problem: 


Maximize  f^xx,  x2,  ...»  xnj 

Subject  to  g^xj,  x2,  ....  xflj  =  b^,  j  =  1,  2,  ....  m;  m  <  n 
Lagrange  method: 

Form  the  function: 


where  the  \-  are  constants  (known  as  Lanra»<ge  multipliers)  as  yet  to  be  deter- 

J 

mined  in  value.  Note  that  when  the  constraints  are  satisfied,  h  is  formed 
merely  by  adding  multipliers  of  zeros  to  f.  Now  treat  x^,  i  =  1 ,  2,  ....  n,  as 
independent  variables,  and  write  down  the  conditions 


3h 


=  0 


0 


ah 

3xn 


=  0 
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EXAMPLE  2: 

Maximize  f(A,B)  =  6A  +  2B  +  AB  -  A2  -  2B2  +  5 

subject  to  p(A,B)  =  2A  -  B  =  8 

Solve  by  finding  the  (local)  optimum  of 


h(A,B,A)  =  f(A,B)  +  A[p(A,B)  -  8] 

=  6A  +  2B  +  AB  -  A2  -  2B2  +  5  +  A  (2A  -  B  -  8) 
Set  partial  derivatives  to  zero. 

3h 


Solving  yields 


Thus 


JR 


=  6  +  B-2A  +  2a  =  0 


|{j-=2  +  A-4B-A  =  0 


3h 

3A 


=  2A-B-8=0 


A  = 


33 

7“ 


B  =  }° 


A  =  1 


.03  TO)  = 


f(^-,  ^-1=  16.5714 
Writing  the  constraint  equation  us  2A  -  B  =  5  g*ves 

2!L .  .  x  -  -l 

36  A  1 


One  can  then  see  that  increasing  6  has  the  effect  of  decreasing  h  (and,  there¬ 
fore,  f)  at  the  rate  of  -1  unit  of  f  per  unit  of  6.  Indeed,  regardless  of  what 
h  looks  like,  if  the  constraint  is  written  as  p  =  6,  then  will  always 

equal  -A. 

Examples  one  and  two  are  both  performed  in  the  same  manner  even  though  one  is 
a  maximum  and  the  other  is  a  minimization  problem.  The  manner  in  which  one  would 
differentiate  between  which  has  occurred  is  by  calculating  the  Hessian  matrix. 


i 

[  APPENDIX  II 

I  QUANTO'S  ITERATIVE  PROCEDURE 

Three  Examples 


This  appendix  is  provided  for  the  reader  to  become  acquainted  /ith  the  types 
of  problems  solved  by  QUANTO.  The  three  examples  serve  to  illustrate  the  three 
basic  options  available  to  the  using  organization.  The  type  problems  addressed 
are 

1.  Optimize  n.. 

*  J 

Given:  three  targets 

three  weapon  groups 

2.  Optimize  aircraft  beddown 

Given:  bomber/tanker  mix 
ten  tankers 
seven  bombers 

3.  Optimize  submarine  locations 

In  cases  1  and  3  the  optimal  n. .  is  found  prior  to  the  optimization  of  the 

•  J 

aircraft  beddown  or  submarine  locations. 

EXAMPLE  1 

Suppose  L  =  3,  M  =  3,  i.e.,  there  are  three  targets  and  three  weapon  groups. 
Let  the  number  of  weapons  in  each  group  be  Nj  =  4,  N2  =  3,  N3  =  7  and  suppose 


0.8 

0.7 

0.9 

s..  = 

1J 

0.6 

0.5 

0.7 

0.2 

0.1 

0.3 

This  matrix  represents  the  survival  probabilities  of  target  i  from  one  weapon  in 
group  j,  e.g.,  the  survival  probability  of  target  2  from  a  single  weapon  in 
group  2  is  0.5.  Let  Vlt  V2,  V3  =  10,  5,  2,  respectively.  As  the  first  step,  an 
arbit;ary  allocation  is  formed,  and  suppose  we  choose 
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where,  for  example,  it  was  decided  to  send  three  weapons  from  the  7  in  group  3 
to  target  1.  The  A  matrix  for  this  (n-jj)  is 


3  n. 


-V"  sn  n  s  • 

3=1  J 


3  n2i 

■V»  s2i  n  s2j  23 

3=1  J 

3  n3i 

'V/n  s„  n  s,.  3J 

3  j=1  3J 


3  nii 

-V"  si2  n  stj  3 

3=1  J 
3  n2i 

-v2fn  s22  n  s2j 

J"** 

3  n3i 

•v3to  s32  n  s3j  J 

J”  * 


3  nii 

s.3  n  s,j  3 

3* 

3  n2i 
-V2£n  S23  0  S2j  3 

3** 

-V3*>  S33  j,  S3j"3j 


and  the  numbers  compute  to  be 


0.82 

0.51 

0.08 


The  smaller  A's  associated  with  target  3  indicate  that  weapons  have  been  over 
allocated  there.  Starting  with  column  1,  i.e.,  the  weapons  in  group  one, 


*31  41  a21 


so  that 


&n  = - - - — 

In  (0.6) (0.2) 


This  ad3'ustment  of  weapons  from  target  3  to  target  2  will  increase  the  objective 
function,  and  will  equate  a31  and  a21.  The  new  (n^)  is 
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and 


Si  nee 


and 


and 


S?53e3i««s* «$£*■ 


2.00 

2 

3 

W2  = 

1.88 

1 

2 

_0.12 

0 

’  2_ 

0.51 

0.82 

0.24 

M2  - 

0.24 

0.32 

0.17 

0.24 

0.34 

0.18 

*22  <  *12 


£n 


0.32 


An  = 


07 


■£n  [(0.7)  (0.5)] 


=  0.88 


N,- 


(Ai  j)  3  = 


Continuing  in  this  manner,  the  Lagrange  multipliers  A.  corresponding  to  positive 

I  X» 

n.&  will  converge  to  the  unique  A&  for  each  weapon  group  i.  The  final  Lagrange 
multiplier  matrix  is 


2.00 

2.88 

3 

* 

1.88 

0.12 

2 

0.12 

0.00 

2 

i 

i 

i 

< 

0.37 

0.60 

0.1  i 

i 

i 

0.44 

0.60 

0.31 

| 

0.24 

0.34 

0.18 

! 
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"0.34 

0.55 

0.16 

0.34 

0.46 

0.24 

0.32 

0.45 

0.24 

which  arises  from  the  optimal  allocation 


4  3 

M =  0  0 

0  0 

which  will  give  an  expect'd  target  value  return  of  14.6  out  of  17,  the  largest 
possible.  This  procedure  .'s  easily  programmable.  For  very  large  matrixes  (on 
the  order  of  several  thousand  targets)  there  are  more  efficient  procedures  to 
adjust  the  multipliers  by  examing  the  convergence  rates.  For  smaller  cases  on 
the  order  or  a  few  hundred  targets  and  weapons,  the  method  abovfe  should  not 
involve  excessive  computer  time. 

EXAMPLE  2 

Suppose  there  are  ten  tankers  and  seven  bombers,  with  a  bomber  twice  the 
value  of  a  tanker,  bedded  down  as  follows: 


where  the  S^'s  are  determined  by  SLBM  yield,  reliability,  trajectory,  *nd  air¬ 
craft  vulnerability,  takeoff  profile  and  sequence.  We  start  by  making  an 
initial  .guess  at  the  allocation  of  SLBMs  to  bases  as  follows: 
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Recalling  that 


and 


U 


2 

1 

0 


3 

2 

2 


S' jB  =  “ViB  (£n  Si jb)  J  SijB 


nij 


XijT 

2 

"ViT  (*■ 

sijT)  n 

S  1J 

bijT 

Xij 

s 

XijB  +  X 

ijT 

0.408 

0.652 

0.193 

XijT 

= 

0.300 

0.407 

0.210 

0.116 

0.155 

0.088_ 

0.161 

0.231 

0.101 

XijT 

= 

0.200 

0.264 

0.147 

_0.029 

0.041 

0.022 

J 

"0.569 

0.883 

0.294 

xi  j 

= 

0.500 

0.671 

0.357 

0.145 

0.196 

0.110 

We  operate  on  this  matrix  column  by  column,  choosing  first  the  column  having  the 
largest  difference  in  A's.  In  column  2,  *12-*32  represents  this  largest  differ¬ 
ence.  The  procedure  requires  a  certain  part  of  the  weapons  in  group  2  to  be 
moved  from  target  3  to  target  1,  since  x32  <  x1£.  Since,  however,  there  are  no 
weapons  in  group  2  allocated  to  target  3  by  our  first  (n-jj)  guess,  we  need  to 
look  further.  The  next  largest  difference  is  .\u-a31,  and  we  move  An  weapons  in 
group  1  from  target  3  to  target  1,  where  An  is  a  root  of 


\ 


aaBMaB^gwaRgsaas 


5** 
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( An  \  _  c  -An  *  4.  c  -An  *  r  +An  ^  ^  ^An  . 

g^n;  -  :>ng  ahb  *  snB  Ai lT  "  a3iB  A31B  '  *3iT  A31T 

This  equation  is  solved  by  use  of  the  Newton  successive  approximation  method  of 
root  finding,  where 

“k+1  =  ink  - 

and  in 

(0.8)'An(0.408)+{0.7)_An{0.161)-(0.3)An(0.n6)-(0.2)An(0.029)  =  0 

An  =  0.88 

Since  g'(An)  >  0  for  bu  >  0,  g(An)  has  a  unique  solution.  Therefore,  our  new 
allocation  is 
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p. 


K 


& 


p 

S’. 

I 

2b 


1 


and  the  value  destroyed  is  about  21.3  out  of  the  total  24.  The  bombers  killed 
turn  out  to  be  5.9  out  of  7  and  tankers  9.4  out  of  10.  An  important  question 
is:  can  the  bombers  and  tankers  be  bedded  down  so  as  to  decrease  the  number  of 
kills  to  a  minimum?  As  shown  before,  the  minimum  damage  that  can  be  inflicted 
by  the  SI  BM  attack  occurs  when  the 


n  s..B 

j=1  ljB 


nij 


i  =  1,  ...»  N 


are  equal  and  when 


are  equal. 


n  S.  *y 
j=1  1Jl 

he  above  example. 


nij 


,  i  =  1 ,  ...»  N 


?  "ijB 


n. . 
ij 


are  0.17,  0.13,  0.12  and 


nij 


4  S«T 

are  0.07,  0.05,  0.06.  These  numbers  indicate  that  the  beddown  is  already  a 
good  one,  and  could  only  be  slightly  improved.  Once  the  given  n^  matrix  is  in 
its  final  form,  in  this  case  with  the  A  -  -  matrix  converged  to  a  tolerance  of 

'  J 

0.001,  the  bombers  and  tankers  can  be  shifted  according  to  the  formula 


/  L  "  *  L  nfcA 

av  =i  n  s  .  -  n  s. . 

\j*1  4=1  kj  J 


applied  to  bomber  and  tanker  values  independently.  With  the  given  initial 
bomber  and  tanker  values  and  a  converged  n.j  matrix  {after  17  iterations) 


3.2649 

3.0 

0.0 

0.7351 

0.0 

4.6912 

0.0 

0.0 

2.3088 
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the  products 


L 

n  s  • . 

j=i ij 


n. . 
ij 


can  be  formed  for  each  aircraft.  Since  there  are  three  targets  in  this  example, 
there  are  three  resultant  products.  These  are 


(1) 

0.1655 

0.0674 

(2) 

0.1289 

0.0547 

(3) 

0.1206 

_0.0621_ 

Bombers 

Tankers 

In  computing  the  values  to  shift,  aV,  for  each  aircraft,  the  smallest  product 
(where  a  value  is  present)  is  subtracted  from  the  largest  one  and  the  difference 
is  multiplied  by  the  value  on  the  target  corresponding  to  the  smallest  product. 

This  is  the  portion  of  value  to  be  subtracted  from  the  total  value  corresponding 
to  the  largest  product.  For  example,  for  bombers,  the  difference  in  maximum  and 
minimum  product  values  is 

0.1655  -  0.1206  =  0.0449 

which  when  multiplied  by  the  value  corresponding  to  the  lower  product  gives 

aV  =  0.0449  x  2  =  0.0898 
Thus,  the  new  bomber  value  matrix  becomes 

8.0898 
4.0000 
1.9102 

By  a  similar  process,  the  new  tanker  matrix  becomes  • 

i 

j 

5.05081  j 

3.9492  J 

1.0000  § 

-1  | 

| 

u 


_ _ _ : _ ^ 


63 


AFWL-TR-73-242 


With  these  values,  the  laydown  is  again  optimized.  By  moving  the  appropriate 

n.  • 

aV  values  iteratively  from  lower  to  higher  n  S^-  ^ ,  we  obtain  the  following 
beddown 


The  new  allocation  of  SLBMs  changes  only  slightly  as  follows 


and  similarly  the  number  of  kills  decreases  only  slightly  to  5.9  out  of  7  bombers 
and  9.3  out  of  10  tankers,  for  a  total  value  destroyed  of  21.2.  Note  that 
although  the  SLBM  allocation  and  total  SLBM  strike  effectiveness  changed  very 
little,  the  beddown,  in  comparison,  changed  considerably.  Hence,  within  the 
context  of  random  events  (i.e.,  SLBM  reliabilities,  CEPs,  etc.)  precipitating 
uncertainties  as  to  exact  numbers  of  bombers  killed,  tne  beddown  problem  does  not 
appear  to  lend  itself  to  a  unique  solution.  Also,  note  that  the  bomber  to  tanker 
value  ratio  of  2:1  was  taken  as  fixed  for  all  bases;  however,  this  value  ratio 
can  be  varied  from  base  to  base.  Also,  note  that  by  inserting  zeros  in  the 
appropriate  locations  of  the  S^.  matrixes,  we  can  effectively  enforce  the  con¬ 
straints  of  not  permitting  tankers  or  bombers  to  be  based  at  particular  locations. 

EXAMPLE  3 

This  example  demonstrates  the  capabilities  of  QUANTO  on  a  small  hypothetical 
problem  which  resembles  those  actually  run  with  QUANTO.  The  input  deck  for  this 
problem  is  listed  in  fiqure  14,  in  the  format  described  in  the  user  documentation 
of  QUANTO.  The  discussion  presented  here  will  be  in  the  form  of  a  guide  to 
reading  the  output,  most  of  which  is  shown  in  appendix  III. 

The  first  two  lines  of  output  serve  to  uniquely  identify  the  run  and  give 
basic  problem  data.  Both  the  beddown  and  the  submarine  locations  appear  on  the 
first  page  of  output.  Targets  1  to  4  are  located  at  points  in  Colorado,  North 
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Dakota,  Iowa,  and  Tennessee  and  the  submarine  locations  are  about  100  nm  off 
the  coasts  of  Virginia,  Louisiana,  and  Oregon.  The  targets  are  distinguished 
not  only  by  their  distances  from  the  submarines  but  by  the  number  of  aircraft 
and  the  distances  (in  nm)  to  the  centroid  of  the  aircraft  flyout  pattern  (a 
function  of  the  numbers  of  runways).  Note  that  aircraft  can  take  off  with  small 
intervals  from  bases  with  dual  runways.  The  missile  parameters  and  trajectory 
data  are  given  on  pages  71  and  72  of  appendix  III.  Page  73  lists  the  input 
aircraft  flyout  profile  data  and  pages  74,  75,  and  76  show  the  profile  generated 
by  QUANTO  with  the  aircraft  leveled-off  at  5000  feet.  The  bottoms  of  pages  77 
and  78  show  the  lethal  radii  and  time  of  (overpressure)  shock  arrival  for  the 
aircraft  of  hardness  indicated  on  page  71,  along  with  other  standard  output  from 
the  nuclear  effects  routines. 

A  table  of  lethal  areas  (as  a  function  of  the  distance  Q  of  the  detonation 
from  the  centroid)  is  then  built  as  each  change  of  distance  to  centroid  (DSPT) 
is  encountered  in  the  target  list.  When  the  aircraft  has  not  reached  its 
terminal  altitude  (at  distance  DISMIN),  new  lethal  radii  are  obtained  based  on 
the  actual  aircraft  altitude  after  it  has  traveled  a  distance  (Q  +  DSPT)  from 
the  centroid,  as  shown  for  two  such  values  of  Q  on  pages  79  to  82  and  85  to  88. 
Occasionally,  a  value  of  Q  results  in  a  geometry  of  lethal  area  in  which  the 
boundaries  of  the  overpressure  and  thermal  kill  regions  have  multiple  intersec¬ 
tions.  When  this  occurs,  approximately  two  pages  of  indicators  are  output  to 
enable  a  detailed  study  of  this  geometry.  This  output  may  be  ignored  on  produc¬ 
tion  runs.  Lethal  area  tables  appear  on  pages  83  and  89  of  appendix  III. 

The  lethal  areas  of  SLBMs  are  based  on  assumed  detonation  points  on  each 
target,  dependent  on  where  the  aircraft  are  at  time  of  weapon  arrival.  Dis¬ 
tances  and  missile  flight  times  from  submarines  to  targets  and  samples  of  the 
computed  aircraft  locations  appear  on  pages  84  and  90  to  92.  Note  that  weapons 
are  numbered  sequentially  through  the  salvos  of  each  submarine  location,  but 
only  the  first  SLBM  of  each  submarine's  two  salvos  appears  in  this  output  in 
order  to  reduce  the  quantity  of  printout  while  providing  enough  information  to 
indicate  when  each  SLBM  arrives  on  each  target.  Note  also  that  the  annulus  P^ 
line  is  only  printed  occasionally;  this  is  because  it  is  not  computed  in 
instances  in  which  the  program  knows  beforehand  that  the  circular  P^  will  be  the 
smaller  of  the  two  P^s.  The  computed  S.^.  values  are  listed  on  page  93  where  i 
is  the  target  number  and  j  is  the  weapon  number,  but  one  row  of  output  contains 
only  the  $. .  corresponding  to  SLBMs  from  a  single  submarine  location. 

*  J 
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The  chart  at  the  top  of  page  94  lists  the  base-by-base  kill  resulting  from 
the  input  laydown  (prior  to  optimization).  Since  each  aircraft  is  worth  one 
point  of  value,  the  number  of  aircraft  and  the  value  are  equal.  The  convergence 
to  the  optimal  laydown  (n-jj  values)  produces  the  output  on  pages  94  to  109,  where 
the  long  form  of  the  output  has  been  requested  to  show  each  An  value,  the 
expected  kill  after  each  shift  of  An,  and  the  allocation  and  multipliers  X-. 

*  J 

after  each  convergence  to  the  tolerances  e  =  0.1,  0.01,  0.001,  and  0.0001.  The 
best  nonintegral  laydown  appears  on  page  107,  listed  by  target  and  then  by  sub¬ 
marine  and  salvo.  The  expected  kill  increased  from  14.9765  to  19.7147  (out  of 
a  total  of  50)  aircraft  during  the  convergence.  Target-by-target  kills  and 
weapons  allocated  (if  n- .  are  rounded  off)  appear  on  page  108. 

*  J 

Optimization  of  submarine  positioning  has  been  requested  in  this  example, 
and  occurs  prior  to  the  requested  beddown  optimization.  The  first  submarine 
move  is  indicated  on  page  110  and  the  resultant  positions  of  all  submarines 
appear  on  page  111  following  the  initial  allocation  (prior  to  convergence  again). 
The  convergence  to  the  best  laydown  with  the  submarines  in  their  new  positions 
follows  with  the  resultant  kill  shown  on  page  134.  The  first  submarine  move 
improved  the  kill  (after  convergence)  from  19.7147  to  20.8742.  The  second  sub¬ 
marine  move  is  shown  on  pages  136  and  137.  A  summary  of  the  submarine  optimiza¬ 
tion  appears  below. 


After  Submarine 


Converged 


Move  No. 

Expected  Kill 

Page 

0 

19.7147 

108 

1 

20.8742 

134 

2 

21.2415 

139 

3 

20.8742 

144 

After  the  third  submarine  move,  the  expected  kill  decreased  slightly  so  the 
submarines  are  fixed  (2  at  point  1,  2  at  point  2,  and  1  at  point  3)  and  beddown 
optimization  begins  with  the  first  shift  of  aircraft  on  pages  145  and  146.  Con¬ 
vergence  to  the  optimal  laydown  follows  each  beddown  change  and  is  summarized 
on  the  following  page. 
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After  Value 

Converged 

Shift  No. 

Expected  Kill 

Page 

0 

20.8742 

144 

1 

18.7694 

148 

2 

16.8932 

152 

3 

16.8171 

156 

4 

16.5884 

160 

5 

16.5863 

164 

Since  the  sixth  value  shift  was  to  be  a  shift  of  less  than  0.05  aircraft 
(see  page  165  of  appendix  III),  the  beddown  optimization  was  terminated.  The 
problem  is  terminated  by  integerizing  both  the  laydown  and  the  beddown.  The 
expected  kill  tends  to  decrease  due  to  the  laydown  integerization  and  increase 
due  to  the  beddown  integerization.  The  resultant  expected  kill  following  both 
integerizations  is  16.4428.  The  inteqral  laydown  appears  on  page  166  and  the 
base-by-base  kills  and  integral  beddown  appear  on  page  167. 

In  analyzing  the  output,  a  table  of  distances  from  submarine  locations  to 
targets  is  useful.  These  distances,  to  the  nearest  nautical  mils,  appear  below 
and  on  pages  84  and  90  to  92. 


Submarine 

Location 

1 

2 

3 

Target  1 

1391 

785 

1107 

2 

1241 

1228 

1174 

3 

881 

895 

1489 

4 

580 

589 

1902 

The  optimal  laydown  from  the  initial  positioning  of  submarines  appears  on  pages 
1 07  and  108.  It  is  interesting  to  note  that  the  submarine  at  location  1  allo¬ 
cated  its  missiles  to  target  3,  leaving  target  4  to  the  submarines  at  location 
2.  Target  2  drew  the  most  weapons  even  though  it  could  not  be  hit  as  soon  as 
the  other  targets.  This  was  probably  because  target  two's  15  aircraft  departed 
from  a  single  runway,  and,  therefore,  were  dispersing  from  a  point  5.5  nm,  in 
this  case,  from  the  brake  release  point.  Dual  runways  permit  more  immediate 
dispersal  since  the  aircraft  can  take  off  in  opposite  directions. 
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The  optimization  of  the  submarine  positions  resulted  in  only  one  submarine 
at  location  3,  even  though  this  location  had  the  best  shot  at  target  2. 

During  the  beddown  optimization,  the  proximity  of  target  4  to  the  submarines 
at  locations  1  and  2  made  target  4  unattractive  for  bedding  down  aircraft  (see 
page  167  of  appendix  III).  The  greater  distance  from  the  coast  outweighed  the 
advantages  of  dual  runways  (i.e.,  immediate  dispersal  due  to  aircraft  taking  off 
in  opposite  directions)  making  target  2  the  most  attractive,  although  targets  1 
^nd  3  also  drew  a  substantial  number  of  aircraft.  In  this  case,  many  aircraft 
would  be  left  on  the  base  when  SLBMs  arrive,  as  is  clear  from  page  91  of  appendix 
III  which  shows  a  zero  inner  annulus  radius  when  weapons  arrive,  even  when  target 
2  had  only  15  aircraft.  Consequently,  it  appears  that  too  many  aircraft  are 
present  to  justify  the  assumption  of  uniform  aircraft  distribution,  and  the  kills 
on  target  2  could  be  underestimated  by  QUANTO. 
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FOR  MlfSllF  TYPF  1 

SUM  OF  LOWEST  l AMROAS  UllH  WEAPONS  =  1.7671 

SUM  OF  HIGHEST  LAMHOAS  =  3.4193 

SUB  MOVEO  FROM  LOCATION  3  TO  LOCATION  1 
THIS  IS  MOVE  NOMHF1'  1  OF  A  SUB. 


110 


\ 


*Z.V  •>  •?*.'<*$$  f^ff*#?***  f&CfitiFZJb  ’-'KWM  ■»^Sv'^*^'"  »’ 


I  I 


M 


•  • 

•  • 

•  • 

• 

A  »* 

«•>* 

r»  >* 

IT  *0 

^  (V 

•4 

z  z 

2  2 

2  2 

Z 

o  o 

O  O 

C  O 

o 

a  a 

a  a 

a  a 

a 

<  < 

<  «x 

<  < 

< 

Ui  UJ 

Ui  UJ 

UJ  UJ 

UJ 

J  > 

*  * 

It  l 

» 

w  w 

««*  w 

w 

-*  nj 

(V 

—  OJ 

t\J 

o  o 

o  o 

o  o 

o 

>  > 

>  > 

>  > 

> 

-J  -i 

-J  -J 

-J  -1 

-J 

<  « 

<  « 

<  < 

< 

w  in 

in  in 

m  in 

in 

*  • 

•  • 

*  • 

• 

oj  rvj 

r>  n 

•H  *H 

tv 

2  2 

2  2 

2  2 

2 

O  O 

O  O 

C  O 

O 

*-< 

M  M 

•»>  ►  * 

M 

►-  *- 

►«  ►■ 

>-  *- 

»- 

<  < 

<  < 

4t  •« 

< 

o  o 

o  u 

a  o 

u 

o  o 

o  o 

o  o 

o 

-J  -J 

-J  _J 

-j  -J 

-J 

CO  £T 

<r  x 

x  m 

X 

3  O 

r>  o 

3 

in  C/) 

in  in 

<n  in 

in 

»-  K 

H-  »- 

H  ►- 

»- 

<  < 

<  < 

4t  <r 

<i 

in  m 

in  in 

in  m 

in 

tr  c 

X  X 

X  X 

X 

in  in 

in  in 

in  in 

in. 

tv  <v 

oj  tv 

cv 

to  V)  l/> 
u*  tu  iu 
J 


©COO 

WOO  4 
*-  o  *  m 
<0^0 
u  •  •  • 
O  to  -* 
-J 


k  *-«  t\j  n 
2 


..vw&ftA&fe* £i 


MvILtlPLIFR  MAT3IX 


> 


n 

o 

< 

o- 

n 

o 

c 

o 

• 

• 

• 

in 

*/) 

IT 

44 

•-4 

M 

'J 

X 

3. 

o 

O 

O 

a 

ft 

ft 

c 

c 

o 

b 

»- 

b- 

* 

n 

* 

CO 

in 

o 

n 

IT 

>c 

>C 

* 

cc 

® 

Jt 

cc 

o 

• 

o 

• 

o 

• 

ft 

o 

ft 

3 

ft 

c 

(Vi 

(VJ 

(V 

44 

x 

(VJ 

z 

n 

z 

* 

c 

* 

c 

* 

c 

ft 

ft 

ft 

u. 

n 

u 

II 

u 

ft 

* 

Q 

(VJ 

a 

* 

a 

ui 

UJ 

UI 

-j 

•j 

-J 

.  j 

z 

44 

z 

44 

z 

44 

a 

z 

* 

a 

z 

X 

ft 

z 

X 

n  ®  * 

O  <\>  ® 

<vi  O'  n 

(vi  cvj  in 

Ui 

CD 

-J 

UJ 

Ui 

CD 

UJ 

Ui 

CD 

3 

t 

1.1 

w^e 

-•  ®  (VJ 

(v  ^  in 

(Vi  (V  (V 

X 

© 

Z 

© 

Z3 

X 

o 

non 

®  O  (V* 

on* 

©  o  *-4 

u 

«J 

3 

u 

_j 

3 

y 

1 

•—  <vj  n  %r 

(VN(C» 

n  (\#  cm 

<  (Vi  (VI  o 

z 

X 

Z 

X 

X 

4K 

•  •  • 

•  •  • 

•  •  • 

•  •  • 

z 

> 

z 

> 

2 

»* 

*“* 

M  44 

#4  #■"• 

z 

*4 

z 

*■4 

z 

•M 

&  O'  *-* 

O'  <0  00 

®  in  in 

(VI  (VJ  o 

c 

•-4 

z 

o 

UJ 

o 

o 

Ui 

o 

•4 

x 

O 

Ui 

*-  tvi  n-  *■* 

t-  m  r*»  -< 

»-  (M  r\j  * 

»- 

K 

H- 

UJ  ®  (VJ  tn 

u;  in  ®  in 

ui  (m  in  :o 

UOOH 

X 

< 

u 

X 

x 

u 

y 

O^MT- 

O  O'  O'  -• 

cm  *  (v 

®  (Vi  (Vi  o 

Or 

b- 

UJ 

ft 

K 

Ui 

ft 

Ul 

ft  •  •  • 

a  •  •  • 

a  •  •  • 

ac  •  •  • 

UJ 

-J 

a 

Ui 

Q. 

UI 

_J 

o. 

<  « 

<t 

x 

►- 

UJ 

X 

»- 

Ui 

X 

fr-» 

Ui 

X 

►- 

►4 

c 

UJ 

*4 

o 

UJ 

*4 

o 

UJ 

116 


DELTA  N  IN  COLUMN  4  FROM  ROW  4  TO  MOW 
EXPECTED  VALUE  KILLED  *  ?0.B7?S 


EXPECTED  VALUE  KILLED  =  ?0.H7*0 


<~^~— -  .  . .  ,  i  i  _ _ _ 


.,  rnyw-vj  **+  te&xiik'gSBe&l I  n''!"'iSM<  ir'Yfl^M  jfl'^Tiilln  vSMli<'"liiif iif'Y7^^^f'i‘"*'rr  IMiiiiiili  i 


EXPECTED  VALUE  KILLEO  =  20.87a* 


H^?V  v^  ..v  .__  ^  ^  _ 


< 

a 

Ui 


'^KKAiMK^ 


X 

ft 

►- 

©  ft  M 

O  MX 

ft  ©  M 

ft  o  o 

~  X 

^'OM 

o  4  n 

O  O  (U 

X 

M  4 

coow 

4*-i n 

4  4 

*-  M  M  4 

MNCtT 

OMm«M 

4  M  M  Q 

•  •  • 

•  •  • 

•  •  • 

•  •  • 

UJ 

H 

M*  M< 

#M  ^ 

-4 

l/>  4  © 

O'  ©  X 

ft  *  m 

ft  ft  4 

& 

HN'Olfl 

^RINN 

hOO  ^ 

•m 

mir.  4  c\j 

uj  in  ©  in 

WNOO 

Ui  4  4 

C  M©  If) 

o  in  in  m 

O  MM  O 

<J 

Q  •  •  • 

&  •  •  • 

ft  •  •  • 

ft  <  •  • 

13 

41  <** 

< 

uc  *■•  *"t 

4[  MM 

X 

►- 

»- 

1- 

© 

© 

X 

•M 

M 

— 

4- 

•M 

* 

» 

2 

O 

c 

o 

ft 

ft 

ft 

c 

o 

o 

p“ 

K 

K 

4 

n 

4 

M 

M 

<VJ 

4 

4 

4 

f- 

K 

P*. 

’B 

® 

* 

© 

X! 

© 

© 

• 

o 

• 

o 

• 

ft 

O 

ft 

O 

ft 

o 

M 

M 

M 

n 

X 

4 

X 

X 

X 

o 

M# 

o 

M 

o 

M« 

ft 

ft 

p< 

ft 

u. 

tl 

u. 

II 

u. 

II 

4 

O 

M 

O 

4 

o 

Ui 

UJ 

UJ 

-J 

-J 

-1 

X 

X 

2 

%Z 

ft 

X 

X 

ft 

X 

X 

ft 

X 

x 

Ui 

Z) 

ui 

UJ 

© 

X 

-J 

UJ 

X 

-J 

Ui 

CD 

J 

Uf 

X 

o 

X 

o 

X 

O 

© 

3 

u 

-J 

© 

u 

~J 

— \ 

u 

_l 

X 

< 

X 

X 

X 

> 

X 

> 

2 

> 

X 

X 

M 

X 

*■* 

O 

a 

o 

o 

o 

c 

M 

X 

UJ 

X 

Ui 

M 

X 

ui 

1“ 

>- 

*“ 

K 

H 

< 

< 

u 

4K 

< 

u 

4 

< 

u 

ft 

K 

UJ 

ft 

K 

Ui 

ft 

K 

UI 

Uf 

-J 

QL 

Ui 

-J 

a 

UJ 

a 

H 

Ui 

X 

f~ 

UI 

X 

►» 

UJ 

x 

** 

o 

Ui 

M 

o 

Ui 

o 

Ui 

128 


EXPECTED  VALUE  KILLED  a  1*0.87 4? 


ssjyjx 


■***:*>■*-> 


I!'  s 


\  A. 


1 

I 

I 

£ 

I 

I 

1 


—  X 

U.  u 


Z  — 

0  3* 

U*  3 
CD  «*  Lu 

TO  3 

-  O  «J 

2  < 

/  > 


n  <  u 

a  k  uj 

u:  _j  a 

u:  x 

— '  C  UJ 


T*  3C 

o  ♦ 

or  o 

t\ 

n  x 
r*  c 
—  a 

u.  it 


(V 

< 

r*> 

»  'O 

o  • 

or  o 

rv: 

x 


a 

fu 

Q 

4- 

o 

J\j 

c 

Ui 

UJ 

UJ 

UJ 

-j 

-J 

-J 

_J 

j 

-J 

•»* 

2 

•M 

X 

M 

2 

k 

a 

X 

k 

a 

X 

k 

0: 

X 

k 

u? 

3 

UJ 

2 

u* 

3 

iU 

ffl 

-J 

UJ 

CO 

Zj 

w 

X 

j 

Ui 

3 

X 

O 

3 

X 

o 

3 

X 

c 

3 

mJ 

3 

w 

-J 

3 

u 

J 

3 

o 

Z 

< 

2 

< 

X 

<T 

> 

2 

> 

z 

> 

z 

> 

Z 

*-» 

z 

*■• 

2 

** 

c 

o 

O 

c 

c 

O 

c 

Uj 

UJ 

*•# 

z 

UJ 

•-** 

2 

UJ 

»- 

1- 

►- 

>■* 

o 

* 

< 

o 

< 

< 

o 

< 

<• 

<J 

lu 

a 

w 

a 

►» 

u; 

a 

►- 

Ui 

a 

w 

-J 

& 

UJ 

-J 

a 

Ui 

.1 

a 

X 

*- 

W 

X 

V* 

UJ 

X 

»- 

Ui 

x 

UJ 

c 

u; 

*■»• 

o 

UJ 

o 

w 

z 

Or  % 
Ut  3 
«  -J 

5  9 

3  1/ 

z 

*  £ 
£  sr 


cxPEcrco  vxtuE  killed  *  so.erug 


/ 


<K>  4*$4S*yvt^WSBaS «•»  jra&*e'*v^?V4  ««<8^^l^c9^C9«X^7«Ei^9W^«^^^,fV%«9T^*^HUH4l^rw.  tiXS*s«i  w.or^««,i'^ 


m  ^ 

©  3 

—  (V 

tv  3 

z  z 

z  r 

z  z 

z  z 

o  o 

o  o 

c  c 

o  o 

ft  ft 

ft  a 

a  a 

a  a 

<  < 

<  < 

<  < 

<  < 

UJ  ui 

Ui  UJ 

ft  ft 

ft  ui 

*  * 

«■»  w 

X  3? 

X  » 

—  tv. 

—  (V 

~  (V 

(V  tv 

o  o 

o  o 

c  o 

o  © 

>  > 

>  > 

>  > 

>  > 

-J  -i 

-J  -J 

ft  ft 

ft  ft 

<  < 

<  4 

41  < 

<  « 

in  in 

in  in 

©  in 

©  © 

•  * 

•  • 

•  • 

•  • 

tVl  (V; 

n  m 

•-4  44 

—  (V 

z  z 

z  z 

z  z 

z  z 

o  o 

°  ° 

o  o 

o  c 

n  r>  <4 

*4  ^  4 

tv  tv 

ft  ft 

ft  ft 

ft  ft 

ft  ft 

<  < 

4t  < 

<  < 

<  •* 

o  u 

a  o 

u  u 

©  o 

ft  ft  ft 

ft  ft  ft 

ft  ft 

o  o 

c  c 

o  o 

o  o 

ft  UJ  ft 

ft  Ui  ft 

ft  Ui 

-J  -1 

-J  -J 

ft  ft 

ft  ft 

o  o  o 

o  ©  o 

o  © 

CD 

aaa 

a  x  cr 

a  a 

m 

CD  CO 

CD  flD 

©  © 

©  © 

■«  <X  <1 

<  <4 

<  < 

»-* 

3  3 

3  Z 

Z  3 

z  z 

ft  ft  ft 

ft  ft  ft 

ft  ft 

in  in 

in  in 

©  in 

©  © 

o  c  o 

o  c  o 

c  c 

*"  ^ 

ft  ft 

ft  ft 

ft  ft 

ft  ft  ft 

ft  ft  ft 

ft  ft 

<  <r 

<  < 

<  < 

4-  tv  tv 

*•  rv  tv 

-4  tv 

in  © 

m  © 

©  in 

©  «n 

<r  m 

©  © 

©  © 

©  tz 

ft 

z  3 

zz 

z  z 

z  z 

Ui 

a  © 

©  © 

©  © 

©  if 

coo 

c  o  c 

c  o 

CD 

>  >  > 

>  >  > 

>  > 

X 

tv  tv 

•«* 

tv  (VI 

(VI IV 

-J  _)  ft 

ft  ft  ft 

ft  ft 

Z 

<  «  < 

<  <  < 

<  < 

z 

©  ©  © 

©  ©  © 

©  © 

X 

X  X 

X  X 

X  X 

X  X 

XXX 

XXX 

r  x 

o 

o  o 

o  o 

o  o 

c  o 

o  c  c 

coo 

o  o 

*4 

ft  ft 

ft  ft 

ft  ft 

a  a 

a  a  a 

a  tr  a 

a  a 

ft 

U.  ft 

ft  ft 

ft  ft 

u.  u. 

u.  u.  u. 

u.  ft  ft 

u.  u. 

< 

ft 

in  m 

©  in 

in  in 

©  © 

©  ©  © 

©  ©  © 

©  © 

UJ 

UJ  ft 

ftft 

w  ft 

ft  ft 

ft  uj  uj 

ft  ft  ft 

ft  w 

ft 

ftft 

-j  -j 

ft-i 

ft  ft 

ftft  ft 

ft  ft  ft 

ft  ft 

»M 

M 

►4  44 

M  M 

•N  M 

•■*  *4 

tn  in 

in  © 

if)  if) 

©  © 

©  ©  © 

©  ©  © 

©  © 

in  m 

in  in 

in  in 

©  if) 

©  ©  © 

©  ©  © 

©  © 

»M  *4 

M  M 

►*  m 

►4  «-4 

W*4  W 

►4  *4  »4 

•*  *"4 

X  X 

X  X 

X  X 

X  X 

XXX 

XXX 

X  X 

z 

z 

z 

• 

—  oft 

tv  ©  © 

DOC 

*  tv  r> 

ooec  ftj 

©  ©  ft  rs 

OO  O 

z 

O  4 

©  © 

©  4 

IT  © 

ftO*  © 

*->o  <  © 

•40  0 

o 

o  ** 

©  © 

©  4 

©  © 

ft  O  >T  © 

>-  o  *4  co 

ft  c  o 

w 

©  ft 

©  © 

©  44 

c  tv 

<0*4  © 

<  O  ft  tv 

<  o  o 

ft 

•  • 

•  • 

•  • 

•  • 

u  •  •  • 

o  *  •  • 

o  •  • 

< 

Kftl- 

ft  44  4* 

*-  (Vi «« 

ft 

O  tv  *4 

O  tV  44 

0*4  4* 

u 

Ui 

ft 

ft 

Ui 

ft 

ft 

ft 

o 

c 

o 

© 

o 

ft 

ft 

ft 

a 

a 

ffi 

© 

© 

ft 

< 

< 

< 

< 

z 

z 

z 

< 

ft 

ft 

ft 

♦- 

© 

© 

© 

133 


3^r^^w5*»*s:  *  r,T^v^v>sxi^ 


c  »- 

C*5 

<VJ 

CO 

a 

tv 

UJ  u. 

X 

© 

K 

n 

-J  < 

© 

sC 

4 

in 

-j  a 

<4 

fVJ 

<7* 

M 

►«o 

• 

• 

• 

• 

• 

x  a 

< 

<30 

<\J 

in 

N« 

© 

i\j 

_J  *- 

© 

© 

o 

o 

© 

<  u. 

O 

o 

© 

o 

o 

»-  < 

© 

© 

© 

© 

© 

C  tt 

c 

© 

o 

© 

© 

1-  O 

• 

• 

• 

• 

• 

X 

M 

< 

in 

in 

in 

in 

50 

C  UJ 

r> 

CVJ 

X 

O' 

UJ  o 

X 

© 

K 

-J  -J 

© 

sO 

4 

in 

-J  < 

-4 

r\J 

fVJ 

•-»  > 

• 

• 

• 

• 

X 

X 

fVJ 

in 

.4 

-J  UJ 

© 

© 

© 

© 

«  3 

© 

© 

© 

c 

>-  -J 

© 

© 

o 

c  < 

© 

© 

© 

© 

»-  > 

• 

• 

• 

• 

IT 

in 

—4 

a.  a 
<  > 
cr  »- 
u 
a 


S'  o 
Z  uJ 

o  *- 

X  4 
<i  u 
u.  O 
r  _• 
j 
<r 


4  rv 


rvj 

-4 

f«* 

CO 

• 

© 

f\j 


II 


o 

UJ 

-J 


e 

5 


•-  a 

U.  Ui 

ex 

o  t 

<  r? 

•-  4t 


iv 


-4 


X 


O 

ui 

V)  r- 

-4  U 


*■ 


»-  a 

o  x 

*-  UJ 


f: 

\ 


t 

4  ' 


134 


fSSL  iki  &„v?  V  Ifarm  .:• 


sC 

c 

X 

o 


z 

o 


II  < 

u 

xr  ~  c  • 
Z  CD 

c  a  z 
a  or  c  </> 
<  •  ►- 
uj  (VJ  < 

*» 


^  it  tvj 

-*  z 

U*  O 


o  x  <  o: 

a  I  C  U; 

7  <  om 

JX 

a  w  3 

>  KI  2 

khitc 
y>  uj  a:  Ui 

itkllU.  > 

JXC  o 
O  O  X 
W  Jl  Ui 
t/>  >  «n 

«-•  fa.  Ik  o 

x  o  o  x 

\r 

ct  x  x  cc  ~ 

03  3  3  X 
U.  </)</>  IT  ►— 


136 


'&■&$$ 


/ 


* 

S’ 


i 

; 


•  • 

•  • 

•  • 

M  .■» 

a  a 

rn  ^ 

.-  (V 

(VJ  ^ 

Z 

o 

M 

Z  Z 

z  z 

z  z 

c  c 

o  o 

o  o 

< 

a  a 

a  a 

CL  Q. 

a 

<  < 

<  < 

<  < 

a 

UJ  UJ 

UJ  UJ 

a  a 

X  * 

X  X 

X  X 

'"■  '■' 

w  «■» 

a 

© 

—  <\J 

*■•  cc 

(\J  rvj 

3 

O 

o  o 

c  o 

o  o 

o 

>  > 

>  > 

>  > 

N- 

-J  -J 

3  -J 

-J  — J 

<  < 

<  •* 

© 

in  in 

in  in 

in  in 

© 

o 

♦  • 

•  • 

•  » 

o 

0j  tvj 

*4  4-« 

— •  CVJ 

o 

o 

o 

© 

z  Z 

z  z 

z  z 

o 

o  o 

o  c 

o  o 

n  n  ^ 

M  «-4  ^ 

M 

•M  urn 

•■4  •-* 

a 

• 

f-  N- 

N-  N* 

K  N- 

a 

<  < 

^  < 

<  < 

> 

u  o 

U  <J 

u  o 

N-  K  ►- 

♦-»-►- 

►-  «—  — 

c  o 

o  o 

o  o 

aaa 

aaa 

►- 

_J  -1 

-J  -J 

-J  -J 

ooo 

©  o  © 

© 

© 

o 

aaa 

aaa 

3 

< 

CD  CD 

CC  CD 

CD  CC 

<r  <  « 

<  ^ 

© 

a 

3  3 

3  3 

3  3 

NNN 

-J 

in  in 

in  in 

in  in 

c  c  o 

ooo 

N- 

h  K 

H>  ►- 

N-  H- 

N»>  N- 

NNN 

< 

<r  < 

M  < 

<  < 

**  (Vl  (Vj 

«"*  (VJ  (VJ 

© 

> 

in  m 

in  in 

©  in 

ff 

© 

ID  CD 

IT  CD 

CD  IT 

3 

0: 

3  3 

3  3 

3  3 

© 

a 

at 

in  © 

in  m 

in  in 

oo© 

ooo 

a 

ID 

>>  > 

>>> 

a 

a 

X 

m  r> 

fVM 

M  n 

-»-J  -J 

3  J  -J 

o 

a 

3 

<<< 

<  <  < 

(vj  r>  © 

** 

Z 

©  ©  © 

©  ©  © 

a 

o 

a 

z 

T  5 

r  a 

X  X 

XXX 

XXX 

or 

N- 

o 

O  O 

o  o 

o  c 

o  oo 

ooo 

X 

© 

a  a 

cr  a: 

a  a 

aaa 

aaa 

3 

3 

N- 

a  u. 

a  a 

a  a 

aaa 

aaa 

Z 

X 

< 

© 

IT 

in  m 

oi  cn 

in  in 

©  ©  © 

©  ©  © 

UJ 

UJ  UJ 

a  a 

a  a 

aaa 

aaa 

a 

»- 

w  -J 

-j  -j 

— j  -J 

jjj 

-J  *j  —J 

a 

M 

M 

M  «M 

•-* 

^4  MMM 

n 

M 

• 

in  in 

©  in 

in  m 

in  ©  © 

©  ©  © 

a 

in  in 

in  in 

in  in 

©  ©  © 

©  ©  © 

a 

a 

►4  M 

M  M 

MMM 

#M  *-4  *"4 

© 

2  2 

2  X 

2  2 

222 

2  2  2 

2 

h- 

Z 

Z 

z 

3 

a 

• 

•-ON 

(VJ 

no  c 

(Vi  n 

o  ©®  (vi 

O©  1*-  co 

o 

Z 

3 

2 

O  4 

O  4 

in  m 

**  ©^  © 

*-©<»© 

NHfyn 

2 

O 

o  — 

©  4 

IT  ft) 

>-  ©*  in 

>-  ©  —  © 

N- 

a 

►* 

O  f-» 

©  — 

©  (VJ 

<  O*^® 

<ON(VI 

« 

z 

H* 

•  • 

•  • 

•  • 

u  «  •  • 

t>  •  •  • 

u 

o 

< 

»-n^ 

V 

NN*- 

h 

O 

cn««^4 

o 

o 

z 

U 

Ul 

UJ 

a 

a 

3 

_J 

-J 

a 

< 

o 

© 

o 

© 

© 

a 

J 

cc 

ft 

a 

a 

CD 

IS 

0 

© 

© 

_J 

< 

< 

■< 

< 

3 

3 

3 

3 

< 

< 

»- 

►- 

*- 

© 

© 

© 

© 

•j 

137  | 


/ 


I 


•  • 

•  •  • 

• 

n  <j 

—  <Vl 

* 

2  2 

2  2  2 

2 

C  C 

o  o  © 

O 

a  a 

a.  a  a 

a 

<i  < 

<  <  < 

< 

UJ  uJ 

UJ  UJ  U' 

UJ 

*  ^ 

XX  3 

x 

w  W 

w  w 

(VI  tvj 

<\j 

o  © 

o  o  o 

o 

>  > 

>  >  > 

> 

-J  -J 

-J  -J  -J 

-J 

«I  < 

<  <  < 

<1 

©  IT 

tr  ©  © 

tn 

•  • 

•  •  • 

* 

cvi  r\j 

t\J 

fu 

2  2 

2  2  2 

2 

c  o 

O  O  © 

O 

r>  r> 

m 

h-  *- 

K*-  K 

*- 

<1  C 

«T  <C  «S 

< 

L>  U 

uu  u 

u 

►-  ►— 

c  c 

o  ©  o 

o 

UJ  UJ 

UJ  UJ  UJ  Uj 

-i  —I 

J  J  J 

-J 

©  © 

©o  ©  © 

tn 

a  x 

xx  x  a 

© 

<C  X 

x  ic  or 

tr 

<  < 

<  <  <  < 

©  © 

©  ©  3 

© 

►—  >- 

in  i/i 

cn  m  in 

tr 

c  o 

cc  o  c 

K  ►- 

*->-*- 

K  ♦- 

►*»  >-  H  ►- 

<  < 

<  *  < 

< 

—  fVJ 

^  fVJ  (V  OJ 

tr  © 

©  ©  © 

tn 

Qr 

rr  2 

©  CO  CD 

w 

u 

tr.  tr 

©  ©  © 

tr 

©  c 

CC  c  c 

i> 

>  > 

>>  >  > 

a 

n  r> 

ftj  wr> 

r> 

_J  _J 

JJJ  -J 

zr 

<  < 

<<  <  « 

2 

tn  tn 

tn©  tr  tn 

z 

J  X 

XXX 

s 

T  X 

a  5  T  X 

o 

c  © 

o  o  c 

o 

c  o 

oo  o  o 

—4 

a  z 

x  *  a: 

X 

X  X 

XXX  X 

»- 

w  u. 

U.U.U 

u 

u  u 

U.U.  u.  u 

< 

a 

©  © 

©  ©  © 

tn 

j.  tn 

©  ©  ©  © 

U.' 

UJ  w 

uj  uj  u: 

u; 

U.  UJ 

UJUJUi  u- 

♦- 

-J  -J 

-J  -J  -J 

^  — J 

-l-l-J  -J 

•— 

M  X 

M  •— 4  M 

** 

•*  M  M 

nj  •><•** «  m 

r» 

• 

IT.  © 

©  ©  © 

tn 

tr  © 

©© ©  © 

© tr 

tn  tn  tn 

tn 

©  tn 

©©  tr  © 

►**  •-* 

M 

** 

W  M 

MMM«4 

X  X 

xxx 

X 

X  X 

_  XX  X  X 

2 

2 

2 

• 

«■*©-♦ 

fvj 

n  ©  ©  <v» 

©  ©  © 

©  ©  4  t\,  4 

O 

2 

©  c 

©  ©  © 

n 

*•©  © 

oo  ©  n 

O 

©  a) 

o  o  tn 

© 

♦-  ©  © 

►“  ox©  © 

*- 

©  4 

c  ©  f\j 

tv 

<  ©  © 

<  o*  ru  tv 

< 

»- 

•  • 

•  •  • 

• 

o  •  • 

U  t  4  •  • 

u 

«« 

k  n 

W 

hM«J 

©  rvj  aj 

o  <*>«-*  *«• 

o 

u 

UJ 

UJ 

UJ 

UJ 

-J 

<J 

-1 

© 

© 

© 

o 

© 

-1 

a 

a 

a 

a 

m 

© 

X 

-J 

< 

« 

3 

3 

*- 

•- 

K» 

in 

© 

© 

138 


vtt,  *>  *.r*c*'r  >/*>*; 


' 


1 

i 


140 


t  -*?y.  ,*j 45*'*^V!rtS5-  •.  &<  fc^*'V  -  vVs^.OTt»v  v 


g*. 


•  • 

•  • 

•  *  • 

• 

.*•  M 

a 

n  <r 

IP  c 

«—  (VJ  *j- 

z  z 

z  z 

z  z  z 

z 

c  c 

o  o 

o  o  o 

o 

a  a 

a  a 

a  a  a 

a 

<  < 

<  < 

<  <  < 

< 

Ui  Ui 

UJ  Ui 

ui  uj  ui 

UJ 

X  X 

X  X 

XXX 

X 

**"’  *“ 

w  '■' 

***  w  *■" 

C\J 

^  Pi 

—  cvj  rvj 

Pi 

o  o 

O  o 

o  o  o 

o 

>  > 

>  > 

>  >  > 

> 

ml  -1 

-J  ml 

jjj 

-J 

<  < 

*«  < 

<  <i  < 

< 

©  to 

to  to 

to  to  to 

© 

•  • 

•  • 

»  ♦  • 

• 

rvj  pj 

ro  n 

-*  PJ 

(Vi 

2  Z 

z  z 

z  z  z 

z 

o  o 

o  o 

coo 

o 

*-*  M 

M  M 

'H  *-»  M 

►-  ►- 

»-  h- 

K 

<■  < 

«  < 

<  <  ■< 

< 

o  o 

o  o 

o  o  o 

o 

o  o 

o  o 

o  o  o 

o 

ml  -I 

-J  -J 

-1  -J  -J 

3 

CL  2) 

CD  X 

x  cc  x 

X 

3  3 

3  3 

3  3  3 

3 

to  to 

to  to 

to  to  to 

to 

h-  ►“- 

t~  *- 

*-*-*-“ 

H- 

<  < 

c  < 

<  <  « 

< 

ip  to 

to  to 

to  11  to 

X 

tt  x 

cd  x 

XXX 

X 

3  “ 

r  3 

3  2  3 

3 

1/  to 

tr  to 

to  to  X 

X 

pj  pj 

MVjM 

(Vi 

UI 

• 

a 

> 

K  H* 

^ 

Ui  UJ 

Ui  UJ  UI  UJ 

Uj  ui 

►- 

©  © 

©  ©  ©  © 

©  © 

X 

© 

a  a 

a  a  a  a 

or  oc 

3 

< 

<*  < 

<  <  <  < 

<  < 

© 

UJ 

H-  ►- 

>-  >- 

ml 

c  c 

o  o  o  o 

o  o 

H* 

K  *- 

< 

(VI 

«-»  (Vi  (Vj  (Vi 

«  (Vl 

© 

>* 

X 

X 

3 

© 

a 

o  o 

o  o  o  o 

c  o 

Ui 

>  > 

>  >  >  > 

>  > 

U. 

U. 

J  J  J  J 

ml  ml 

© 

a 

<  < 

<<<< 

<  < 

(Vi  (VI  -4 

*-« 

to  © 

©  ©  ©  © 

©  © 

a 

o 

UJ 

3  X 

X  X  X  X 

I  X 

X 

K 

o  o 

o  o  o  o 

o  o 

X 

© 

*  a 

x  a  a  a 

X  X 

3 

3 

u.  u. 

a  a  u  u. 

a  u. 

Z 

X 

©  © 

©  ©  ©  © 

©  © 

© 

uj  u; 

UJ  Ui  UJ  U) 

Ui  UJ 

X 

ml  ml 

ml  ml  ml  mJ 

3  mi 

Ui 

M  M 

lyi  M  M  N 

X  •-.  *-» 

►4 

1  ( 


no  ci\j 
c  c  -c 

©  O  IP 
O  OftJ 


IP  ©  ©  10 

<p  IP  ©  © 

M  M  M  H 

X  X  X  X 
X 

OO^ftJ4 
e  o  o  n 
H-  O  «0  IP  O 
<  O  •*  (VJ  0J 
<_)•••« 
O  (V  ~ 


XX  X 

Z  3 

o  ©  o  z 

MOO 

►-  o  ©  »-  -•  (M  n 

*  ©  ©  z 

u  •  •  »-» 

o  ^  o 

-J  CL 


_  » 

u&y&v  A  1  i&£  ±*£k 1  ■* 


Vv  ■*■ 


nBMnRi 


■ffrcs*  *v  *  - 


(\J 

• 

z 

o 

M 

H- 

< 

a 

Ui 

h- 

>?  run  ** 

►< 

o  ^  n  ® 

UJ 

*  «D«  #- 

t/) 

•  •  •  •  u» 

3 

Ul 

•*C 

o 

o 

Q 

o 

►- 

►«* 

MD  II 

C 

o 

o 

* 

e 

o 

US  US  US  Ui 

o 

o 

o 

a  a  a  a 

o 

>>>►*- 

o 

*••  Uj 

«*• 

O 

m 

►-►-*-»-  ct 

u.  u.  u.  u.  * 

^  *  <  <  ♦- 

ana  i 

f- 

O  O  U  U  X 

tfl 

a  a  a  a  o 

■« 

•—  •-  •-  •—  a 

Ui 

<  <  <r  <  U. 

-J 

•  •  •  •  c 

► 

— •  n.  r-  ^  Uu 

c 

> 

c 

V 

X 

X 

K  K  ►  H  #• 

a 

fc!  Uj  U  Ui  •" 

Ui 

O  <£  O  O 

u. 

a  t  to: 

u 

<  <  ^  < 

•-  «-  •“  »**  Ui 

c 

a  r 

t  o  a  o  >  u- 

*- 

C  C  C  C  *-  X 

tr. 

U.  U.  u.  U.  X 

3 

U.  3 

T 

O  Z 

O  <J  o  o 

i/> 

::2:kh 

a 

CCCCkk 

u* 

C  C  C  O  41  ~ 

aa  3  aa  2 

_j 

a  a  a  a  c  j* 

a 

a 

+m 

>~  V  V  >-  —  U.’ 

•- 

^  <  3 

-J 

-J 

2? 

JJ  J  Jlf  < 

X 

> 

-a  t  a  rr  o 

ik 

<  <  <  <  it  «n 

o 

>  >  >  >  •  •- 

z 

•■»  — •  *-»  ++  r*> 

«* 

>  >  >  >  «r 

a 

a  a  a  a  *-• 

o 

=  333  x 

< 

<r  <r.  it  v>  »- 

-j 

145 


i 


*3^^i5X>e^S^SK  .^feSS^*'-* O**  S*0'“"^vA^v<^V 'W!  ;-r 


OK 

<\j 

4) 

4) 

o» 

Uj  U 

« 

ac 

K 

n 

K 

-J  < 

o 

CO 

* 

K 

ja 

* 

K 

<Vl 

4 

•"•  O 

• 

• 

• 

• 

• 

x  tr 

nn 

< 

<r 

<\J 

4fl 

*-« 

CO 

-J  K 

o 

o 

m 

o 

<  U 

o 

<4 

o 

i/) 

o 

K  - 

o 

o 

o 

a* 

o 

O  . 

o 

IT 

o 

o 

K  U 

• 

• 

• 

• 

• 

a 

tn 

X 

•-4 

15 

05 

c  u: 

r* 

uj  3 

X 

-J  -j 

o 

-j  < 

—•  > 

• 

X 

<\J 

CC 

CO 

K 

r> 

CD 

■fr 

K 

h» 

O' 

(M 

• 

• 

• 

f\J 

in 

«-* 

< 

a 

UJ 


-J  UJ 

C 

If. 

o 

IT 

< 

c 

* 

c 

XT 

K  _j 

c 

© 

o 

c> 

o  < 

c 

IT 

© 

K  > 

• 

• 

• 

• 

IT 

a 

if 

o 

o 


K  UI 

u  a 
<  > 
a  k 
u 
a 


*/*  o 

Z  uj 
O  H» 
X  < 

<  o 

uj  O 
"X  -J 
-J 
< 


K  O' 

u  UJ 

o  » 
a  x 
«  z 


H*  z 


146 


z 

o 


< 

a 


UJ 

<u 

irn  c  ir 

^c  ft  —  x 

-4NX1 

UJ 

4  T>  IT  4 

tr 

•  •  •  «H 

2 

u. 

« 

o 

U 

a 

< 

2 

H- 

►- 

H  il  it  ll 

C 

c 

^  — «  #-«  — «  »-> 

o 

e 

*-* 

o 

o 

o 

uj  u.  u.  u; 

c 

a  a  a  a 

e 

>>v>k 

c 

H  K  f»  h  l^j 

_ 

c 

• 

a 

u.  u.  u.  u.  < 

<<<<*- 

a  or  a  or 

h- 

U  O  U  U  X 

x  ft  a  or  c 

< 

a 

U- 

<  <  <  <  u. 

-J 

•  ♦  ♦  ♦  c 

h- 

«-  ft  r.  4  u. 

< 

> 

o 

> 

X 

X 

H  fr~  »-  K  ~> 

* 

U.  U.  lu  u.  rv. 

u. 

o  c  o  c 

u. 

Iftft  1 

u. 

<  <  <  < 

M 

^  ^  t“  V-  Uj 

c 

a  a: 

ft  a  a  a  >•  x 

o  c  O  C  K  £C 

tr 

u-liLU.  T 

2 

U.  2 

X 

*-  »-*-*-  C  Z 

u  u  o  o 

(f) 

3  2  2  2  *«  K 

ft 

C  C  C  C  U.  It. 

u. 

C  O  O  C  <  •-• 

a  x  a  a  a  x 

-J 

a  a  a  a  o  tr 

ft 

ft 

>-  >  >  >  ~  U' 

►— 

hi-  <  r 

-J 

Jj J JC< 

X 

—  —  fh  > 

a  oc  x  cr  a 

u. 

<  <  <  <  4  Ji 

C 

>  >  >  >  •« 

z 

—  >o 

< 

>  >  >  >  ^ 

ft 

ft  ft  ft  ft  •-« 

e 

2  2  2  2  X 

< 

tf  tr  tr  tr  t- 

-j 

149 


• 

•  •44 

•  • 

•  • 

mm 

M  M  M  A 

a  >~- 

<<•»  a 

o 

n  «*  cp  *.o 

«  <VJ 

<v  * 

2 

ZZZ2 

2  2 

Z  2 

O 

o  c  o  o 

O  O 

o  o 

a 

a  a  a  a 

a  a 

a  a 

< 

<  <  <  < 

<  <. 

<  < 

UJ 

UJ  UJ  UJ  UJ 

UJ  UJ 

UiUJ 

j  m 

X  X 

X  * 

** 

■*»  ^  *■■* 

W  V 

**  w 

^4  (V.-iOJ 

•—  CVi 

rvj  rvj 

o 

o  c  o  o 

o  o 

o  o 

> 

>  >  >  > 

>  > 

>  > 

ml 

•  -J  -J  -J 

•J-i 

-J-l 

41 

<4  «  «  «C 

41  4| 

41  < 

cp 

to  cr>  cn  cr» 

<p  cp 

CP  CP 

• 

»  •  »  • 

•  • 

♦  ♦ 

PVJ 

r\j  <v  tn  m 

N  —4 

-4  C\J 

2 

2222 

2  2 

2  2 

c 

c  c  o  o 

O  O 

C  O 

r> 

*-  cvi  (\i  -X 

(V;  (Vi 

M 

mm  mm  mm  mm 

mm 

•—  mm 

►- 

»-»-*-  H- 

H  H- 

h”  K* 

4  4  4  4 

<  4T 

•4  < 

o 

uuuu 

u  a 

a  o 

i- 

h- 

o 

o  o  o  o 

o  o 

o  o 

UJ  UJUJ 

UJ  UJ  UJ  UJ 

UJ  UJ 

-1 

J  J  J  J 

-j  -j 

— J  -J 

o  oo 

o  o  o  o 

o  o 

*— 

a  acr 

a  a  a:  oc 

Or  QC 

c 

CD 

tr  cc  ma 

X  X 

CD  CL 

<  <  « 

4  4  4  < 

<  < 

4- 

3 

3  3  33 

3  3 

3  3 

»- 

H*  •- 

CP 

CP  CP  CP  CP 

CP  CP 

CP  CP 

o  oc 

c  o  o  o 

c  c 

►- 

K  H- 

*-  *~ 

KKk 

K  -  KK 

< 

<  <  «J  •* 

<  41 

4T  < 

(UCVJ 

*4  -4  (VJ  CVi 

pm  (V. 

CP 

CP  CP  CP  CP 

CP  CP 

CP  CP 

a 

CD  CD  CC  CD 

Cf  X 

CC  CL 

a 

* 

3  3  2  3 

3  2 

3  3 

u 

CP 

CP  (P  CP  CP 

CP  CP 

CP  CP 

C  C  3 

c  c  c  o 

c  c 

® 

>  >> 

>>>> 

>  > 

(VI 

(Vl  (V.  — *  »— 

CVi  (V 

(V  (Vl 

-J  «J  -J 

J  J  J  J 

-J  ml 

“■* 

4  44 

4444 

<  * 

2 

CP  tP  IP 

CP  CP  CP  CP 

CP  CP 

- 

2 

X  X  2  X 

X  X 

X  X 

XXX 

X  X  X  X 

X  2 

C 

O 

c  o  o  o 

c  o 

o  o 

coo 

c  o  o  c 

c  o 

•M 

OC 

ir  ec  a  cc 

a:  x 

oc  a 

a:  croc 

oc  oc  a  oc 

a  cr 

•>. 

U. 

kkiklk 

u.  u 

u.  u. 

u.  u.  u. 

u.  u.  u.  u. 

U_  U- 

2 

CP 

(P  cP  tp  CP 

IP  CP 

CP  Ip 

CP  CP  CP 

IP  CP  CP  CP 

CP  CP 

u 

w 

k.  UittJUi 

x  k 

UJ  U. 

UJ  U*  UJ 

W  UJ  UJ  u- 

UJ  u * 

»- 

-J 

_J  _J  _J  _J 

-J  -J 

-J  -J 

-J  JJ 

_J  -J  -J  _l 

-J  -J 

pm 

MMMM 

mrn  —• 

M  *4 

mm  mm  pm  m 

ru  *-  mm  mm  mm 

r.-w 

CP 

CP  CP  CP  CP 

CP  CP 

CP  CP 

CP  cp  to 

CP  CP  <p  cn 

CP  IP 

CP 

cr  cr  cp  *p 

CP  CP 

CP  CP 

CP  CP  CP 

tr  cp  ip  cp 

CP  CP 

MtMMM 

mm  mm 

M 

mm  mmmm 

mm  Pm  mm  mm 

►4  pm 

X 

X  X  X  X 

X  X 

X  X 

XXX 

X  X  X  X 

T  X 

2 

2 

2 

• 

«—  e 

rv  o  c  ©  © 

r  or 

4NO 

c  o  nt* 

c  o  o  c  c 

C  c  o 

2 

r* 

ref  oo 

C  4 

IT  IT 

O  *CP 

•-r-  r>  ip  m 

e  c 

O 

«— * 

cc  cc  o  c 

c  ip 

KOlf-4 

—  C  CO  — 

HCC 

4V 

necoo 

©  4) 

D  — 

<  o  4n 

4  ^  n  cc  •■* 

400 

*- 

• 

•  •  •  • 

4  • 

•  • 

O*** 

O  •  •  •  4 

U  •  • 

< 

•»  — 

h-  pm  mm  mm 

o  oj 

04-  - 

C  mm  mm 

u 

Hi 

UJ 

UJ 

UJ 

J 

JJ 

ml 

c 

O 

O 

o 

O 

-J 

a 

oc 

a 

ar 

CD 

OD 

CL 

-1 

4t 

< 

4f 

41 

3 

3 

3 

< 

»- 

f- 

*- 

CP 

CP 

CP 

151 


,5^1 


-  -  •"* 


• 

<f  •  •  • 

•  • 

• 

«««  M  ^*> 

«« 

(•3 

n  4  in  x 

—  rvj 

(VJ 

2 

z;:zz 

z  z 

Z 

O 

oo  oo 

o  o 

o 

a 

a  a  a  a 

a  a 

a 

< 

<  <  <  < 

<  < 

< 

Li 

LJUJUJ  UJ 

UJ  Li 

UJ 

X 

X*  *  * 

X  9 

V 

w*  w  w  w 

w  M* 

w 

—  «vi  — ■  <vi 

—  (VJ 

(VJ 

o 

OO  O  O 

o  o 

o 

> 

>  >  >  > 

>  > 

> 

-J 

J  J  J  J 

-j  -j 

-J 

< 

<  4  «4  4 

<  < 

«x 

c n 

x  x  x,  x 

C/>  (/V 

x 

» 

»  •  »  » 

•  » 

» 

(M 

(VMVr  tn  r> 

M  M 

m 

z 

z  z  z  z 

z  z 

z 

c 

oooo 

O  o 

o 

r>  m 

4 

^  (V  (Vj 

CVJ  CVi 

M 

M  (M4 

M 

>• 

H-  H* 

H-  H* 

H- 

< 

0404 

<  < 

< 

o 

oooo 

o  o 

o 

H-  h- 

►- 

K  K  K 

h  H- 

o 

oooo 

o  o 

o 

Ui  UJ 

UJ 

UJ  UJ  UJ 

Li  UJ 

_J 

J  J  J  J 

_J  «J 

-J 

o  o 

© 

©  O  © 

©  © 

03 

a  ct 

X 

x  x  a. 

X  X 

O 

tt 

a  o  x  Qj 

X  X 

X 

<  < 

4 

4  <  4 

4  4 

4 

o 

oooo 

o  o 

o 

►-  H* 

K 

H-  *- 

1/3 

x  x  x  x 

X  X 

IT 

c  c 

o 

o  o  c 

c  o 

H- 

HKKH 

*- 

K 

K  K 

K 

f-  *“ 

K  H- 

•c 

<4  4  4 

«4  4 

« 

m*  rvj 

(Vi 

M*  M 

(V) 

x 

x  x  x  x 

x  x 

x 

X 

XX  X  X 

X  X 

a 

a 

oooo 

o  o 

z 

Li 

x 

X  J/3  IT  1/3 

in  x 

in 

c  c 

c 

COO 

o  o 

03 

>  > 

> 

'  >  >  > 

>  > 

7 

(Vi 

(vj  ru  — •  -« 

CVJ  (Vi 

(Vi 

— j  -J 

-J 

-J  -J  -J 

-j  -J 

o 

4  « 

4 

444 

«  4 

Z 

in  in 

X 

XXX 

XX 

r 

X 

T  J3  3f 

X  X 

z 

J  z 

X 

X  X  z 

X  2 

c 

o 

C  O  O  O 

o  o 

c 

oo 

c 

o  o  o 

o  o 

a 

x  x  a  x 

a:  x 

X 

X  X 

X 

XXX 

X  X 

►- 

L 

L  U.L.  L 

L  L 

L 

L.  L 

L 

L.  L  U» 

L.  L. 

4 

or 

x 

x  x  x  x 

in  in 

in 

in  in 

X 

XXX 

XX 

Li 

L> 

u;  Li  u;  uj 

UJ  UJ 

Li 

Li  Li 

ui 

UJ  Li  Li 

UJ  Li 

1— 

-i 

.j  j  j  j 

«J  -j 

•j 

-J  -J 

J  J  J 

••  _J 

••* 

M  MM  M 

M  M 

#M  «M  Mt 

M 

IV  M#  «M  M* 

r. 

X 

x  x  x  x 

x  in 

in 

mir 

X 

XXX 

X  X 

X 

x  in  in  x 

X  X 

X 

iiii/1 

X 

XXX 

XX 

•m 

MMMM 

M  M 

m 

M*  M 

•M 

M  M  M 

M*  M 

2 

Z  Z  Z  Z 

z  z 

z 

z  z 

z 

X  Z  X 

z  z 

z 

z 

z 

► 

**  <3 

(V;  4  o  ©  a 

n  o  (Vi 

4  X 

O  o  (Vi 

X 

C  ©  4  O 

o  O  O 

2 

O 

OOOO 

O  -« 

<© 

tM  O  «M 

OD 

IM  O  X  o 

MOO 

O 

o 

oooo 

O  (M 

r* 

w  O  <Vi 

(*•* 

H  X  o  O 

►-  o  o 

HI 

1/3 

4  0  0  0 

o  r*. 

(Vi 

<ON 

(VJ 

4  X  4  o 

4  OC 

w 

• 

•  •  •  « 

•  • 

■ 

O  *  e 

• 

u  •  •  • 

o  •  • 

< 

K*  (VJ  *— *  • 

►-  (VJ  -* 

K 

o  (M 

O  *"•  (M 

U 

LU 

UJ 

UJ 

Li 

-1 

•J 

J 

o 

c 

© 

O 

«j 

x 

X 

X 

a 

X 

03 

QD 

.j 

< 

< 

< 

c 

o 

O 

o 

*4 

*. . 

*- 

X 

X 

X 

155 


** *«WK*iw,i»i(wi  titr&gG;*£ns9B  &2K#*iamr»* 


w 
u.  a 
<  > 
a  ►- 
u 
a 


i/>  c 

2T  UJ 
C  W 
£.  < 
<  CJ 
UJ  c 
V  -J 


CO 

'O 


o 

Ui 


»-  a 

UJ  UJ 

O  X 
a  3 
«  ID 


< 

»- 

O 


UJ 

3 

-J 

> 

O 

UI 

*- 

o 

UI 

a 

x 

UJ 


156 


1*5*60  MJSSILFS  FPOM  ?  SI/HS  AT  SUP  LOCATION  2.  SALVO  1  (WEAPON 


N 


(V.  D  **  |fl  <0 


2  2  2  2  2 
o  c  o  o  o 
a  a  a  a  a 

<  <  <  <  < 
Miumnui 
*  *  *  *  * 


f\j  ~  (\j  r\j 


o  o  o  o  o 
>>>>> 
-J  -J  -j  -j  -j 

<c  <  <  < 

in  in  in  in  in 


•  • 

• 

—  (V 

(V 

2  2 

2 

o  o 

o 

a  a 

a. 

<  < 

UJ  UJ 

LU 

9  * 

* 

—  f\i 

tV 

o  o 

o 

>  > 

> 

-J  -J 

^  < 

m  m 

in 

iv  iv  m  m  — 


2  2222 

2  2 

2 

c  o  o  o  o 

C  O 

O 

r<  a;  «**,  ^ 

^  (Vl  (V. 

tv.  tv 

W-4 

*-* 

♦—  *- 

>— 

<  <  <  <  < 

<  < 

UUUL'U 

o  u 

V 

H  Kl-K 

K  K  K 

»•  »— 

c  c  c  c  o 

c  o 

© 

UJ  UJ  UJ  Uj 

UJ  UJ  UJ 

UJ  UJ 

-J  -J  -l-J  -J 

-J  -J 

- J 

O  O  O  O 

©  ©  C 

©  o 

o:  a:  a  c  <r 

CD  X 

a:  a  a  a 

a  a  a 

a  a 

cc 

<r  <  -fl  <r 

■*  «  < 

<x  <r 

O  O  o  33  O 

O  O 

D 

KKK 

h  k 

IT  (T<  t/)  IT 

m  m 

in 

C  C  C  C 

COO 

C  O 

*-  H* 

»- 

K  KKl- 

»-  h 

<  <  <  <  < 

it  in  m  it  i/i 

*c  < 

< 

—  tv  tv.  tv 

«f«*  (Vj 

~*tv 

in  it 

in 

x  i  x  ec  a* 

x  a: 

a? 

O  15 

=> 

IT  IT  IT  1/1  m 

in  i/i 

/> 

o  c  o  c 

c  c  o 

c  c 

<Vi  tv  fv:  — 

>  >> > 

>  >  > 

>  > 

IV  CM 

(V 

J  JJJ 

-J  -J  -J 

— J  -J 

<  <r  <  < 

«  <  < 

<  < 

it  in  in  in 

ip  in  sn 

m  in 

s  r  :  j  y 

5  a 

O  o  c  o  c 

o  c 

a  a  a  a  a 

a  a 

u  u  u.  u  a 

U-  lv 

iSi  tn  ts>  m  m 

m  m 

Uj  U  u»  uj  u* 

U.  UJ 

j  J  ^  j  j 

->  -J 

in  m  in  m  in 

in  m 

in  <r.  in  in  in 

:r  m. 

M  M  M 

M  M 

HUT 

X  X 

tv  in  o  c  o  o 

nocr 

€C  ^  O  O  O 

o  a 

jitoco 

o 

(\i-tooo 

o  ^ 

•  •  •  •  4  •  • 

h  «  (Vi  —  ►-  <v 


u 

UJ 

UJ 

UJ 

c 

c 

e 

o 

_J 

a 

a 

a 

-J 

< 

< 

4K 

< 

►- 

P- 

3 

2 

3  3  X 

2X3 

3  3 

C 

C 

o  o  o 

o  c  o 

oo 

a 

a 

a  a  a 

a  a  a 

a  a 

u 

u 

u  u  u 

u  u  u 

u  u 

in 

»n 

vi  m  tn 

'/)  v'i  n 

in  vi 

w 

U 

Uj  U  UJ 

U*  UJ  UJ 

U*  UJ 

mJ 

-1 

-J  -J  «j 

~J  .  J  —J 

-J  -J 

*•* 

^ 

tn  *«4  »-• 

m 

m 

ip  m  m 

m  in  m 

<n  m 

in 

in. 

in.  in  in 

in  m  m 

m  in 

«»NM 

M  M 

«  *«• 

X 

X 

XXX 

XXX 

X  X 

2 

2 

2 

^  p- 

c  c 

in  a  p* 

C  o  o  c 

C  e  e 

►“  c 

a  a  — 

-«'JC 

o  e 

a 

p-  c 

a  ~cc 

P^iTC 

>-  •  o  o 

IV 

<  o 

tV*3  IV 

<IP*JC 

<  cc 

« 

©  • 

•  •  • 

(J  •  •  • 

o  •  • 

K 

otv 

«■* 

o  —•  tv 

c 

UJ 

-J 

_J 

J 

o 

a 

X 

m 

a 

© 

o 

o 

m 

in 

m 

159 


TARGET  WEAPONS  AIRCRAFT  TOTAL  KILLED  TOTMu  KILLED 

NOMPEP  ALLOC* TEO  TYPE  VALUE  VALUE  AIRCRAFT  AIRCRAFT 


o 

o> 

O' 

C\J 

tr 

a 

« 

o 

o 

® 

O' 

IT 

O' 

e 

• 

• 

• 

• 

a 

cv 

PVJ 

— 

<\ 


(V. 


o 


-a- 

a 

a 

IT 

>C 


r\.  m  4- 


cr> 

-J 

<r 


c 

w 


c 

Ui 

»- 

o 

Ui 

a 

x 

Ui 


160 


4 

V 


. .  •i'*»«*S 


z 

o 


< 


a 

UJ 

*«• 

*- 

a:  rvj  <•  « 

M 

ft*  f-  * 

f"»  <0  « cu 

UJ 

X  <£ 

If) 

•  •••♦- 

r» 

UJ 

<x 

O 

u 

c 

< 

o 

*~ 

II  II  II  II 

c 

o 

«  — •  "*  !• 

c 

c 

ft 

o 

o 

c. 

U.  U.  U.  U. 

o 

a  a  a  & 

o 

>►►►♦- 

c 

1“  »"  h  K  !i; 

,_, 

O 

• 

*-*-►-►-  a 
u.  u.  u.  u.  < 

<  <  <  <  K 

i  >- 

U  U  u  U  T  l/i 

i  saao  < 

*-  —  *-—  ^  uj 

<  <  <  <  ,*.  _j 


•  •  •  *  c 

—  rv; <-  ^  u. 

c 

> 

c 

> 

■*■ 

or 

►*  K  I-  H  — 

a 

U  u.  U  u  if 

u. 

oocc 

u 

i  a  3-fi 

u 

<  <  c  <r 

aj 

c 

a.  a 

o  .•>  >  u 

*— 

c  c  O  w  *-  T 

tn 

u.  U.  U.  u.  T 

u  z 

T 

^  ©  2 

uuuu 

•J) 

r  z  -  r  ►■  ►• 

'X 

c  o  O  C  U  U 

UJ 

c ccc<— 

a  o  a  a  a  x 

-J 

a  :  a  a  it 

a 

a 

v  >  >■  y-  +•  ix. 

*- 

h»  •-  h-  ^  2 

-j 

*“•  ♦»  M  1 

D 

JJ  JJN<* 

X 

—  *•*  ♦-  •-  C  > 

T  t  ff  I  «C 

Uj 

<<  <<  Oifl 

© 

>  >  >  >  *  •-• 

X 

*»■«  *—•  M  M 

< 

>  >  >  >  l/> 

ft 

a  a:  ft  ft  ^ 

o* 

zzzz  X 

< 

*/*  IT  IT.  if  *- 

-J 

161 


»  aajassaaagggiflaaaBaaaaiBaiaaBQagssitti  gasattaaaaaaa5e3BM«isaaiaffifiaiiaBesaaaai^aaa^MBBaB^«aat£aa5aiaaa5 


£S2SSSsS32S325SK^ 


^  <v  c  c 

IT  vC  -« 
fw  ^  ^  ^ 

•  •••>- 


u:  u.  u,  u. 
a  a.  a  a 
>>>•>■»- 
»-*-►->-  u. 

o 

•-►•Kwfl: 
u.  U.  U.  Ll  « 
«  <  <  ^ 
a  a  a  a 
U  U  U  U  2. 

or  a:  a:  arc 

•“•**.»■•«•“<  Q 
*»  <  <  <  U. 

•  •  •  •  c 

— •  (V  f*>  ^  w 

> 

c 


U  U  CM 
3  c  ^ : 
o  c  c  ( 
c  o  3  t 
a  a  a  < 
a  a  a  t 


-  err 

x 

.  -  c 

>  X  ui 
fsi 

-  ►*  «■*  c 
.  u.  a  u.* 

l  »"  UJ  N 

:  2  C*- 
MTU;a 
»-  Ui 

•  us  T  o 

t  3  ~ui 


J  j  jjit<22 

M  MMM  ^  >  O  •*"* 

tt  AC  CD  QD  «■*  •— 

<"  <  <<Ci/)»*2 

>>>>  •  tM  <  } 

u  o 

>  >  > >  woe 

e  a  ft  a  h jo 

ZJki 

in  in  ir  m  h  <  tr 


(\j  V  tP  <£  <n  (Vj 


222Z 

c  o  o  c 
a  &  a  a 
<  <  <  < 
Ul  UJ  UJ  UJ 


M  VMM 

•«p  — <■ 

*M 

(\)  (VjmM 

—  rvi 

o 

o  o  o  o 

oo 

> 

>  >  >  > 

>  > 

-J 

— I  _J  — J  -J 

-J  -J 

< 

<  <<  < 

<  « 

CP 

lp  tn  tP  IP 

IP  ip 

» 

•  •  •  * 

»  » 

(VJ 

^  cvj«n  n 

#M  *-* 

2 

2X22 

2  2 

o 

o  o  o  o 

O  O 

n  cvj  «p 

«-»  (Vj 

f\J  (VJ 

M 

M  M  M  M 

M  M 

»-  c- 

< 

«t  <  ^  < 

<  < 

u 

U  U  U  U 

u  u 

»-  ►- 

»-  »- 

»-  »- 

o 

o  c  o  o 

o  o 

UJ  UJ  UJ 

UJ  UJ 

UJ  UJ 

-J 

-J  -J-J  -J 

-J  — 

o  o  o 

o  o 

©  o 

O' 

X 

a  ICCDD 

a  or  a 

a  a 

a:  cr 

(*> 

DC  QD 

<  <  < 

<«  41 

<  c 

D 

3  3  3  3 

3  3 

»-  >- 

K  ►- 

IP 

IP  IP  IP  IP 

IP  CP 

coo 

o  o 

o  o 

K 

t-  H- 

*-  »- 

K  H 

< 

<t  <  <r  < 

«  < 

~  <v»  cm 

(Vj 

«  fVJ 

ip 

CP  IP  CP  (P 

if  iP 

a,  CD  CD  IT 

j-t  if.  cp  cp 


J  2  I  J 
C  O  C  C 
a  a  tr  or 
u.  u.  u  u 

IP  IP  IP  (P 
UJ  UJ  UJ  US 
J  JJ  J 

MM  M 

ip  cp  €/>  tr 
cp  it  ip  tr 


o  o  o 
>  >  > 
-J  -J  -J 
<  <  < 

<p  tr  tr 

o  o  o 
or  ac  a 
u.  w  u. 

<p  ip  tp 
UJ  UJ  U.' 
-J  -J  -J 

I  M  M  M 

to  IP  IP 
cf  CP  cp 


x 

X 

XXX 

X  X 

XXX 

X  X 

X  X 

z 

2 

2 

o 

(VJ  c 

oce 

POC 

* 

o  o  o  o 

o  o  o 

coo 

o 

c 

©GO 

o  e 

►-CO© 

MOO 

M  o  O 

o 

o 

©  O  © 

©  o 

*-  ©  ©  © 

>-  ©  © 

►-  ©  © 

c 

o 

©  ©  C 

©  © 

<600 

<oo 

<oo 

• 

• 

•  •  • 

•  • 

u  •  •  • 

CJ  •  • 

a  •  • 

h-  ev 

K  *“• 

(VJ  mm  ^ 

H*  (Vj  m 

i- 

o  (VI  ^  ^ 

O  Pj  (Vj 

O  MM 

UJ 

UJ 

UJ 

tal 

-J 

-J 

_J 

c 

© 

© 

O 

a 

a 

a 

ac 

'  0 

X 

QD 

< 

< 

< 

« 

3 

3 

3 

*- 

*- 

CP 

CP 

CP 

'CtoUoMMES u. 


AFWL-TR-73-242 


REFERENCES 


1.  Eisen,  Dennis,  and  Alderson,  Linda,  "Advanced  Bomber  Pre-Launch  Survivability, 
Lamda  Corporation,  Arlington,  Virginia,  June  1972. 

2.  Everett,  Hugh,  III,  "Generalized  Lagrange  Multiplier  Method  for  Solving 
Problems  of  Optimum  Allocation  of  Resources,"  Operations  Research,  V.  11, 
pp.  399-417,  1963. 

3.  Hiller,  F.  S.  and  Lieberman,  G.  J.,  Introduction  to  Operations  Research, 
Holden-Day,  Inc.,  September  1968. 

4.  Thomas,  M.  A.  and  Gemmill,  G.  W.,  Increasing  Salvo  Kill  Probability  Through 
Aim  Point  Patterning,  US  Naval  Weapons  Laboratory,  Technical  Report  (NWL) 
TR-2'756,  bahlgren,  Virginia,  June  1972. 


